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COMPUTERS IN THE AEROSPACE CORPORATION 



I . Introduction 

In the modern aerospace corporation, computers are used widely in the 
scientific, engineering and other technical operations. Probably more so than in 
any other industry, they have been accepted as key contributors in the equipment 
research, design and development activities and beyond these they play increasing 
roles in manufacturing, test and check-out, and field operations. Computers are 
now so commonplace in the typical company that in contrast to the situation 15 years 
ago when one of the pioneers of the industry thought that a few computers spread 
around the country would satisfy all possible needs, one wonders whether a few com- 
puters spread around a company will satisfy the needs within that company. With 
them playing a dominant role in the technical activities it is useful to review 
their functions and clarify by reason of these functions, the qualities to be sought 
within the computer and the manner of its use in research, design, and manufacturing. 
Rather than having a pre-disposition to this computer or that computer we shall look 
at general purpose computation and by relating the desired result with the possible 
computer mechanization derive a look ahead to the future aerospace computer depart- 
ment . 

Technological breaks- through in the components of computers are not ex- 
pected within the next decade. Solid state circuitry is commonplace and wherever 
it provides the desirable performance it is presently in use. High speed memory 
devices in digital computers and widebandwidth computing components in analog com- 
puters, equal to the tasks they perform are available in every manufacturers equip- 
ment. The improvements to be anticipated are in man-machine communication and com- 
puter organization, the first to improve the ease and flexibility with which the 
computer is applied, the second to increase computing capability and value. Both 
improvements require detailed attention to be paid to the way computers are used 
in solving useful problems, that is to say we must look to where, why and in what 
general manner rather than to the detail of individual programming. Fortunately, 
there is a consistency of methods and procedures in the many aerospace companies 
that allows us to do this. We can categorize and classify, and this procedure is 
very useful and informative. 

II . Aerospace Systems Design and Development 

Looking at the progressive development of a system in an aerospace com- 
pany, be it a shoulder- launched , rocket-propelled, infra-red homing missile for use 
against low-flying aircraft, a super-sonic variable wing, intercontinental airliner, 
or a manned orbiting space station intended for one year's successful operation, 
the ste-by-step procedures are similar. Perhaps slight variations enter and possibly 
there are different useful arrangements of the list but the eight steps to success- 
ful operation of the system are : 

1. Conceptual Studies aimed at sketching out the overall performance to 
be attempted, the state-of-the-art components, the organization and integration of 
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the system, the major difficulties to be attacked and solved, the alternative ideas 
that can usefully be investigated when all of the way-out possibilities have been 
legitimately discarded. Conceptual studies are performed by men with ideas, re- 
searchers and designers both, but those very much interested in the future rather 
than the present. They produce typically a paper proposal of a system never be- 
fore built, not able to be built today, but certainly by projecting today's labor- 
atory experimental device into tomorrow's mass-produced item, something practical 
in the future. Until one begins to build it, the limitations of pencil and paper 
must be supplanted by the almost unlimited boundaries of the computer to assist in 
the testing of ideas. 

2. Experiments to establish basic physical data. Typically a projected 
aerospace system is touching the bounds of the unknown at several points. The 
strength of alloy materials at high temperatures, the effects of fatigue on plastic 
materials, the structural integrity of a thin variable plan-form wing under gust 
loading, the sloshing in a liquid-fuel rocket and its effect on vibration and thrust 

control, the aerodynamic forces at high Mach number, the expansion of a turbine 
blade at high speed and temperature, the spectrum of visible and invisible light 
that is useful for guidance, the effects of cross-coupling non-linearities on atti- 
tude control, etc. etc.. Basic experiments must be performed to establish the data 
for use in the design of a total system and these experiment's both in their planning, 
their efficient execution, and in the appropriate reduction of the measured data to 
a more usable form benefit from the use of computers. 

3. Component Design . A total system is the integration of a set of sub- 
systems made of components. From the projected performance of the system one places 
requirements on the behavior of the sub-systems and thereby the components. A first 
step in the practical realization of the system must therefore be the design of com- 
ponents; gyros, accelerometers , hysteresis motors, clutches, synchros, gear trains 
for the homing system of a missile; pumps, turbines, tachos, variable ducts, combus- 
tion chambers, fuel injectors, for the jet engines of a supersonic transport; thermo- 
couples, humidity devices, gas analyzers, air conditions, air purifiers and all of 
the other components of an environmental control system for a manned orbiting space 
station. In some instances, component design may not be necessary for many com- 
ponents will be available from other systems but wherever it is necessary, the 
modern computer has a part to play in establishing desirable parameter values for 
optimum performance according to recognized criteria. 

4. Sub-system Performance Analysis . The assembly of components into a 
subsystem (a propulsion system, an attitude control system, a guidance system, a 
structural system, etc.) always introduces the possibility of poor interaction 
between the components. This is particular true where energy feedback loops exist, 
permitting dynamic instabilities and completely unacceptable behavior. The indivi- 
dual components can be treated as "black boxes" having transfer relationships be- 
tween input stimuli and output signals, as established in the component design phase 
of the program, and their assembly into a working sub-system is thereby amenable to 
overall performance analysis, usually making use of a computer. 

5. System Integration and Hardware Check-out . As the system design moves 
along and pieces start falling into place in the form of real hardware, the dynamic 
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behavior of the total system, made up as it is from the many sub-systems, comes in- 
to focus. The test of overall performance before all parts are available is a stan- 
dard procedure to avoid a "back- to- the-drawing-board" situation when the first pro- 
totype system "flies". Computers, of course, play key roles here for theoretical 
analysis and "seta-of -the-pants" performance conjecture is quite inadequate. It is 
at this point that the concept of the engineering simulator comes through in big 
terms, permitting a real-time representation of parts of the system not yet avail- 
able in real hardware to be used in the exercise of existing hardware. The homing 
missile control sub-system hardware can be exercised in a realistic target approach 
situation, the transition phase in the variable wing plan-form can be exercised be- 
fore flight under control of a simulated auto-pilog, the attitude control system of 
the manned orbiting space station can be tested just as though it were in space and 
without the assistance of aerodynamic damping. Beyond this engineering simulator 
role, computers assist in the solution of many problems in the organization of the 
total system and its eventual production. 

6. Mission Planning . Aerospace systems, be they for research, commer- 
cial transportation, or military purposes are always intended to move from one point 
to another. The missile is launched and approaches the target, the supersonic trans- 
port flies from New York to London, the orbiting space station is placed in orbit 
and then must be re-supplied by satellite vehicles. There is typically a require- 
ment, therefore, to plan these missions so that by some criterion of performance 
they are completed successfully -- the missile must destroy the target no matter 
from which point it is launched, the transport must consume minimum fuel and main- 
tain its schedule through a bad-weather landing condition, the supply vehicles 

must rendez-vous successfully with good energy management. In all cases the point 
of interest has now moved from the design of the system to its effective exploita- 
tion in its projected environment. Analyses of this kind are, to be sure, beyond 
simple tasks of calculation and are now possible only through the extensive use of 
computers. Quite naturally, one of the most important missions to be planned is 
that of the test flight, or normally the many test flights, to ensure effective 
instrumentation and flight coverage in the proving of the system. 

7. Crew Training . Commonplace in the aerospace industry is the con- 
cept of training simulators, by which computer models of the designed system and 
the environment in which the system operates are used to train the staff who are 

to operate the system. For the all-important reasons of cost and time, crew safety, 
and instructional convenience, training simulators have become a significant part 
of the overall aerospace systems' business. Once they used special purpose com- 
puters, but now they depend almost exclusively on general purpose eouipment frequently 
equal or similar to that used in the engineering simulator for system integration. 
The infantryman must fire simulated missiles, the pilot must fly a simulated air- 
plane, and the crew must control a simulated space station because for training 
purposes it is cheaper, more efficient, safer and convenient. 

8. System Check-out, Flight Testing and Evaluation . The final step in 
the design of any aerospace system is taken by a demonstration of satisfactory per- 
formance. The many ranges in this country and others attest to the frequency with 
which aerospace systems are checked out tested. In every case the large volume of 
data that flows from tracking and telemetry stations presents a formidable task in 
computer data reduction and evaluation. Without modern general purpose computers 
the task could not be completed. 
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These eight phases of the modern aerospace system design process require 
computers at every step. The tasks they perform are very varied and it is clear 
that the demands placed on the computer, the mathematical operations performed, are 
peculiar to each task. Despite this, the methods and procedures of modern design 
cause such tasks to be separable into three groups and these groups result compu- 
tationally in a need for (a) Simple Calculation, (b) Data Reduction and Signal 
Analysis, or (c) System Simulation. 

(a) Simple Calculation . The tasks of scientific research and engineering 
design have been continually increasing in complexity. However, many times a task 
requires relatively simple calculations involving algebraic operations and perhaps 

a few transcondental functions. When the volume of work warrants it, this drudgery 
is now transf errable to small calculators, readily programmed, and rapid in their 
operation. They range in size and cost from PDS-1020 computer ($20K) to the IBM 
7094 ($3M). 

One would expect to find the need for simple calculations throughout the 
program reviewed above, but particularly in phases 2 and 3. 

(b) Data Reduction and Signal Analysis , as the name implies, recuires 
the computer to accept input data on which it acts to provide a refinement, either 
a reduction in sheer volume with the establishment of important parameters that 
are equally valuable (mean, standard deviation, correlation coefficients, power 
spectral density, equivalent system constants, etc.) or a tranformat ion to a form 
more useful to the engineer who needs to use the information, or yet again the es- 
tablishment within the data of some key signal (a natural frequency, a mode shape, 
a sequence of events, etc.). The task is typified by relatively simple mathemati- 
cal processing of fairly large volumes of data. It suggests relatively unsophis- 
ticated computation but a need for appreciable speedy operation and good machine 
communications. In conventional hardware, the computer could be either analog or 
digital but because data reduction and signal processing has typically been an off- 
line operation with a consequent need to store the information, it is now widely 
performed with digital computers the signals to be processed being sampled and 
stored on magnetic tape. It is interesting to speculate on the value of a hybrid 
computer for such tasks for this has both the data storage capability of digital 
equipment and the signal filtering qualities of the continuous analog devices. 

One might expect to have tasks in data reduction and signal analysis in 
phases 2, 3, 4, and 8. 

(c) Systems Simulation . The major capability provided by modern com- 
puters in the activities of aerospace systems design, the one that has changed the 
procedures of design over the last twenty years, is that of system simulation. The 
ability to model on a computer the important behavior of any physical system and 
manipulate this behavior at will by adjusting first the system parameters and then, 
if not satisfied, the structure of the system until acceptable behavior is obtained 
is extremely valuable. Before any item is produced, it is thereby possible to 
"destroy a target", "fly to London", "orbit the earth", at least, the ideas of how 
to destroy a target, how to fly to London, how to orbit the earth can be tested. 
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For system simulation purposes, the computer requirements contrast with 
those of data reduction. Here the need is for relatively sophisitica ted mathemati- 
cal programs acting on very small amounts of input data (system parameters, condi- 
tions, and environmental stimuli). Conventionally the task has been one matched 
by the speed and program sophistication of the parallel analog computer. Analog 
simulators have "flown" many aircraft and missiles into combat before the prototype 
left the runway. More recently the increasing speed of digital computers, though 
painfully slow in comparison, has caused simulation programs to be developed for 
the more capable of these machines. 

One might expect to apply simulation procedures in all phases of a sys- 
tems design with the possible exception of experiments to establish basic physical 
data. The capacity required from phase to phase varies and the methods of exploiting 
simulation change but computer simulation is so universal that it deserves consid- 
erable discussion both to place it in perspective and to derive the consequent com- 
puter characteristics. 

Ill . Design by Computer Simulation 

Question: What is the sequence of the design process at any point of the 
development of an aerospace system? 

The answer to this question, at least for a high proportion ot situations, 
is well demonstrated by a simple block diagram of a form familiar to control sys- 
tems engineers. The design process begins with a statement of reouirements , usually 
incomplete, sometimes supplied by another agency, department, or section, which es- 
tablishes a required performance. We can look upon the statement as a block in the 
diagram and the signal at the output of this block can be viewed as the required 
performance, P r . 

The signal P r is applied to the engineer who acts in the form of a compar- 
ator with error-detecting qualities. He applies -- a corrective signal, E, which 
one can look upon as an excitation, to a block which we can describe as creative 
design. This block, powered by two supplies, knowledge (of the useful physical 
phenomena) and experience (of the past successful and unsuccessful uses of these 
physical phonomena) synthesizes a design which can be represented as the informa- 
tion or signal, D, at the output of the block diagram. 

Having produced a design, the process works in a feedback manner to test 
the design. The feedback loop contains a block which we can title analysis for 
with a proposed design it is possible to analyze performance, Pq;, in order to check 
design performance against required performance, P r . Any difference results in an 
excitation, E, stimulating the creative design. Hopefully we have a stable, though 
high-gain, system and obtain successful results rapidly. Occasionally the process 
is unable to bring E to zero and adjustment is necessary in required performance, P r . 

In this analogy between the design process and a feedback control system 
the content of the block termed analysis is informative, if not a little amusing. 
Three kinds of feedback element are possible : 
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(1) Mathematical analysis, characterized by an inexpensive two step pro- 
cess of (a) problem description and (b) theoretical evaluation. The problem descrip- 
tion is today a wide bandwidth, unity gain, filter able to pass adequately the quali- 
ties of a complicated and possibly sophisticated system. The theoretical evaluation, 
however, is a narrow bandwidth, high order filter unable to handle any high frequency 
qualities of system sophistication. The feedback element, therefore, provides a 
rapid but severely limited evaluation of performance, obscuring the sophisticated 
qualities . 

(2) Physical experimentation, characterized by a typically expensive wide 
bandwidth filter in cascade with an extensive transport delay, passes the qualities 
of a system for evaluation very well but only after considerable time. This element 
demands the very low gain in the creative design forward path that existed previous 
to World War II, for stability to be maintained. The result is very slow progress. 

(3) Computer simulation provides an optimum feedback element for it places 
a short circuit between the more desirable parts of mathematical analysis and physi- 
cal experimentation by-passing the theoretical evaluation and the transport delay of 
building experimental equipment. It represents an inexpensive two-step process of 
mathematical model building and "computer" experimentation which provides for a 
rapid, high gain design process even when complicated systems are to be designed. 

With this introduction simulation as an essential part of modern aerospace 
design and as the most significant use for the modern computer, we face the question 
of what kind of computer matches the needs of today and tomorrow. Perhaps one com- 
puter designed to fulfill the simulation function will also qualify as a data pro- 
cessor and calculator. Perhaps not. 

IV. Computers for Simulation 

A computer which is to be used for simulation purposes must provide facility 
for building a structure which corresponds in space and time (or at least gives the 
impression outwardly of corresponding in space and time) to the system and environment 
simulated. It must have the added qualities of permitting a simple building process, 
a simple injection of desirable stimuli and the effective measurement of consequent 
behavior, and for the future it must permit the development of automatic design pro- 
cedures in order to reduce time-consuming "cut-and- try" methods of today. 

If one considers the typical use of a design simulator one recognizes three 
parts to its programming which are essential to the experimental procedure. Firstly, 
there is a physical system to be represented. This is the subject of design and its 
representation must allow adjustment of its structure and parameters. Secondly, to 
understand and thereby evaluate the behavior (performance) of the system when sim- 
ulated it must be placed within an environment which is also programmed into the sim- 
ulator. Thirdly, the evaluation process must be controllable by the injection of 
desirable stimuli with consequent measurement and the possible use of automatic de- 
sign up-dating. 

Although individual examples have particular program demands for these 
three parts, sometimes containing little environmental program and at other times 
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an extensive one, an example is useful to illustrate this program division. Con- 
sider a simulation aimed at designing an automatic landing system for the super- 
sonic transport. It has features similar to those found in the simulation of the 
homing maneuver for the guided missile, of a space vehicle rendezvous and docking 
maneuver, and many other situations in which two bodies are recuired to meet in 
space with an appropriate attitude but no specification of time. From the point 
of view of the simulation program what do we find in the three parts? 

The physical system probably contains relative position and velocity 
sensing devices (a radar possibly), an autopilot for attitude control, a jet en- 
gine thrust controller, and the various aerodynamic control surfaces with their 
local position controllers. The environment, in this example, contains the aero- 
dynamic forces and their consequent integration through the equations of motion 
to establish aircraft position and attitude. The stimuli are both predictable 
(changing initial situations as might result from consideration of the local air 
traffic) and random (wind gusts, radar noise, etc.). The measurements of perfor- 
mance are landing force and velocity, and the accelerations experienced at the 
pilot's seat during let-down. 

This example suggests the kind of programming that is needed in aero- 
space simulation, and the qualities (amplitude and frequency) in a real-time re- 
presentation , 



For many years now analog computers have been used with great success 
to simulate physical systems, especially those which were the subject of engin- 
eering design. It is likely that they will continue to be used with equal success 
in the future. They are characterized by favorable qualities for this kind of 
work, the most important of which should be stated in a discussion of possible 
computers for simulation. 

(a) Analog computers are readily programmed for simulation by the 
parallel interconnection of elementary computing devices to match the ecuations 
that describe the original system. 

(b) Analog computers obtain a relatively high speed computing capability 
from their "one device per required operation" organization which permits a con- 
current computation of many parts and thereby typically a "faster than real time" 
simulation of common aerospace systems and situations. 

(c) Analog computers contain an extremely simple method for integra- 
tion with respect to time, facilitating ready solution of nonlinear differential 
equations and thereby a simple means for simulating the time behavior of physical 

systems . 



(d) Analog computers represent the variables describing a system to be 
simulated by continuously variable voltages which are then readily displayed in 
parallel on commonly-used strip chart recorders, oscilloscopes, and x-y plotters. 
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(e) Analog computers have a "get-at-abili ty" which is considered by most 
scientists and engineers to be desirable when employing an equipment for simulation 
purposes. One has a feeling of contact with the system being represented, through 
the "one-to-one" relationship between sub-systems and computer hardware and also 
through the ease with which program rearrangements and system parameter adjustments 
are obtained. 

(f) Analog computers havea resolution equal to or better than the needs 
of most engineering systems to be simulated. However, care must be exercised in 
parts of problems which are of increasing importance in the aerospace industry. 

Does this list of favorable qualities omit any that are desirable in com- 
puters for simulation? In other words in the years during which analog computers 
have been used in simulating systems like the automatic landing system would any 
additional ability have been useful, or can we foresee an increasing number of sit- 
uations where other qualities would be desirable. 

Consider the automatic landing problem and we can draw some conclusions. 
Firstly the generation of aerodynamic derivatives of one, two, or possibly three 
variables is inconveniently accommodated by the common analog function generators. 
A difficult problem to solve in any manner but we can seek alternative methods of 
information storage. Secondly, the possible need to incorporate aircraft structural 
behavior due to its interaction with the attitude control system suggests at least 
a more economical method of representing distributed systems. In the same vein, 
the completely parallel nature of the analog computer organization may not be solely 
a desirable feature for it implies an unreasonable large number of components being 
required for complicated simulations, a number which can rapidly get out-of -control . 
Thus we might seek to incorporate in the computer a method of using components se- 
quentially at high speed, similarly to the use of the arithmetic unit of a digital 
computer. 

Next we are almost certainly today in the simulation of complex aerospace 
systems to find need for representing sampled-data devices, and logical elements. 
More and more digital devices are being incorporated in such systems and thus a 
simple computer mechanization of their characteristics is necessary. Finally, the 
resolution with which variables of primary importance to an investigation can be 
represented is increasingly insufficiently handled on the analog computer. It is 
true that rescaling is possible, even automatic rescaling, but it would be much 
more acceptable if the computer could represent an approach and touch-down maneuver 
from 10 miles range and 2,000 feet altitude without any changes. 

We can therefore suggest adding the following abilities to the computer 
for simulation. 

(a) An ability to store data in a flexibly changeable manner so that sys- 
tem parameters and environmental functions can be readily loaded into the simulator 
and adjusted at will, even automatically under the control of the computer. 

(b) An ability to perform computations at high speed and sequentially in 
order to economize on the use of analog components that are assigned to exactly 
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similar operations (e.g. in the representation of a distributed system having impor- 
tant structural, thermal, acoustic, or electromagnetic behavior outwardly the simu- 
lator would appear to represent the physical system faithfully in space and time, 
even though only sample points of the structure are displayed.) 

(c) An ability to simulate logic systems described by discrete-valued 
binary variables rather than continuously-valued variables. 

(d) An ability to provide when necessary an increased resolution in the 
simulation of a situation. (Homing studies, Space-vehicle trajectories, etc.). 

To these abilities in which the conventional analog computer is deficient 
can be added two very desirable computer qualities that are not well-represented in 
the analog computer. 

(e) A convenience in applying the program to the computer, checking it 
out, and storing it for future use, (in place of the present mechanical patch panel, 
etc.). This would reduce the turn-around time for the general purpose computer. 

(f) A self -programming quality, or at least in place of it the possibil- 
ity of "softwre" rather than hardware which reduces the tasks of programming. 

Not unexpectedly, these desired additional abilities and qualities are 
well represented in the conventional general purpose digital computer. Its three 
prime features are (1) flexible program and data storage in random access, high 
speed memory or slower speed bulk storage, (2) an economical sequential computing 
procedure which uses a common operating unit to perform many mathematical operations, 
and (3) a method of computation which provides increased resolution, theoretically 
without limit at the expense of time for computation. The general purpose digital 
computer has been constructed conventionally in a manner which does not provide the 
favorable qualities listed above for the analog computer. It is therefore reason- 
able to attempt to match the requirements of modern aerospace simulation by a deriva- 
tive of the two computers -- a hybrid computer. 

To be universally valuable in the role of a simulator a hybrid computer 
cannot be a simple combination of an analog computer and a digital computer, each 
allowed to communicate with the other. It must be a derivative of the two equip- 
ments (use both continuous voltage and binary word representations of variables 
which describe the physical system), of the two organizations (operate in parallel 
and serially), and of the two computing methods (be programmed so that units can 
be operated sequentially and/or concurrently as best suits the requirement s) . 

To be practical, the hybrid computer of today will take advantage of ex- 
isting hardware. It will be modular in construction, allowing the build-up of 
capacity with time from a basic unit. It can be viewed as having five parts: 

1. An analog section with a parallel array of high speed integrators, 
summers, multipliers and function generators. 

2. A digital section with logic and arithmetic ability, and volatile 
data storage. 
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3. An analog-digital data communication system to translate analog 
voltages to binary numbers and vice-versa. 

4. A program control capability to maintain the timing of events in 
the computer. 

5. Input-output display and eouipment monitoring for program checkout. 

A typical use of such a system would find the analog section performing 
high-speed integration of non-linear differential equations with the digital sec- 
tion providing precise low speed computation of slowly changing variables. These 
variables might require the numerical integration of differential equations or 
might be the result of an algebraic function generation. Any other function gen- 
eration or algebraic transformation might fit on the analog equipment or be per- 
formed digitally. Samples of computed data would be stored digitally and logical 
data evaluation would also be performed digitally. The program timing and control 
would be a responsibility of the digital equipment. However consistent with their 
normal separate operation, an application would require all operations on the dig- 
ital computer to be completed cyclically, whilst the analog computer might have two 
modes of operation, one a high speed iterative mode, the other a continuous compu- 
tation mode. 

The kind of operation envisaged for such a system, does not demand an ex- 
pensive digital computer. As might be expected in the combined use of analog and 
digital computers, the 3-4 decimal place accuracy of the analog computer (~ 12 bits) 
is well matched by a 16-32 bit range in word length on the digital computer. This 
requirement appears to follow through no matter the type of digital section consid- 
ered, be it gpds, DDA, special purpose or multi component digital computer. The 
digital section needs to be fast, for the parallel, high-speed nature of the gen- 
eral purpose analog section is otherwise frecuently wasted. Other features of the 
digital section are the subject of considerable discusiion. 

The choice of linkage equipment is rather easily made. The accuracy of 
the analog section determines the word length required in the converters. Typically 
it is 12-13 bits plus sign. 

A number of A-D channels are required. Signals can be simultaneously 
scanned, stored in analog track-and-store units, and converted sequentially by one 
converter to binary number form for storage in the digital section. Otherwise they 
can be sequentially scanned, without the need for analog storage, and converted to 
binary number form. The D-A channels are either separate, each channel accepting 
binary numbers representing consecutive values of a variable and converting them to 
analog voltages which are updated at the operating rate of the converter, or they 
also can use a multiplexer. If a multiplexer is used then again the analog signals 
are stored between the updating times on track-and-store devices. The conversion 
rates depend on the settling time of the converter. The D-A converter operates in 
a parallel manner and can be designed to reach a steady value in no more than 10 
microseconds. The A-D converter operates serially, the diode gates settling one 
after another, and can be designed to reach a steady state in approximately 60 micro- 
seconds. Thus analog signals can be translated into binary numbers at a rate of 
15,000 samples per second. 
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The control operations desirable in a hybrid computer must cover the needs 
of typical applications. Certainly there is need for being able to control the mode 
of operation of the analog section rapidly and consistently. The operation of in- 
dividual integrators or their operation in groups must be controlled from computed 
signals. The intention to exploit the high-speed capability of the analog section 
implies that all switching operations desirably be electronic. Relay switching but 
presumably will slow the simulation of some systems. The sampling of computed values 
in the analog section must be controllable. Read-out and read-in of data should be 
automatically controllable. This suggests automatic coefficient setting equipment 
under the control of paper-tape, digital memory, or computer values, and a complete 
monitoring system similar to these already available on modern analog computers 
(ADIOS on the EAT. 231-R system). 

V . Mathematical Characteristics of Design Problems Requiring Hybrid Simulation 

From experience with problems to be solved in the design of aerospace sys- 
tems one can suggest that problems having one or more of the following mathematical 
features should lead to a consideration of hybrid computation for system simulation. 

1. Simultaneous differential equations with widely different parameters 
which produce both low and high frequencies in the solution. 

2. Differential equations to be solved at high speed, their solutions 
for different initial conditions or parameters being used in a prediction, itera- 
tion or optimization process. 

3. Combinations of continuous and discrete variables as in the descrip- 
tion of a sample data system or a computer control system. 

4. Perturbation analysis about slowly changing, precisely established 

solutions . 

5. Statistical analysis requiring repeated solution of differential equ- 
ations, including Monte Carlo methods for deterministic problems. This is essen- 
tially a data storage and simple evaluation task around the solution of differential 
equations . 

6. Filtering and processing continuous and sampled data for evaluation 

purposes . 

7. Partial differential equations to be solved by serial integration 
procedures . 

8. Ordinary differential equations accompanied by transport delays. 

A few examples of practical applications of hybrid computation in aero- 
space system design help to clarify the situations in which these mathematical fea- 
tures occur. Consider the problem of designing the guidance and attitude control 
system for a space vehicle attempting to dock with a target vehicle, possibly the 
manned orbiting space station. This docking maneuver might be manually controlled 
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by the pilot in which case the design of the control mechanisms must take due con- 
sideration of the problems of a man-machine system. A close-range maneuver control 
system design suggests immediately the use of an analog computer to simulate the 
intended system. The choice of the analog computer is particularly emphasized by 
the desirability of checking the design with a man "in the loop" in real time. Al- 
so the selection of control system parameter values, and the possible re-arrangement 
of any use of sensors for feedback purposes is particularly easy when they are mech- 
anized on such a computer. The later inclusion of real hardware is more readily 
accomplished using the analog computer. However, there are features of such a con- 
trol system which are discrete or sampled and, therefore, suggest the use of digital 
circuitry. The reaction jets which adjust the attitude of the vehicle are "on-off" 
rather than continuously adjustable. The manual 3-axis control stick used by the 
pilot energises in a number of pre- selectable ways some, but rarely all of, say, 16 
jets. Between the control stick and the jets there is a logic unit whose operation 
is modified by other manual inputs. These two features of the system could be mech- 
anized with difficulty on a conventional analog computer. They are readily established 
with the logic capability of digital equipment. 

The same maneuver control system could be completely automatic with a bea- 
con-homing device replacing the pilot's decisions, but in either case, the system is 
most likely to many parts which are discrete in their operation and thus call for 
digital equipment in their representation. 

A second example of the desirable use of a hybrid computer is in the anal- 
sys of rocket parameters and staging times for maximum efficiency in orbit injection. 
A multi-stage rocket's flight into orbit with the necessary considerations of aero- 
dynamic forces, on-board control system, and staging with its consequent sudden 
changes in mass, inertia and thrust, is a complicated problem in simulation. To 
achieve any efficiency in considering many possible system parameters and flight 
path designs, a high speed simulation of a typical six minute flight to injection 
is essential. Thus, the speed of parallel computation of an analog computer is 
desirable. The automatic control and evaluation of the many different sets of para- 
meters with their effect on the selected criterian of efficiency is most appro- 
priately achieved using a digital computer. Such an arrangement would save con- 
siderably on the time presently required by an analog computer and even more so on 
that of an IBM 7090. This kind of problem, a search for an optimum set of parameters 
for a dynamic physical system, is one having many varied application. Its implications 
in computation - the solution of differential equations plus the need for up-dating 
stored values and making of logical decisions according to a criterian - suggest 
strongly the use of a hybrid computer. 

There are many other examples of problem types which place a mixture of 
demands on computers. They should be most readily solved by the use of a hybrid 
computer. In choosing such a computer there is good reason to consider both the 
technical features essential in the computer, and the economic comparison between 
the different possible organizations. 



Hybrid Seminar 
Page 13 



VI. Hybrid Computer Applications 

Although the different parts of a hybrid computer can be usefully applied 
to simulation tasks by themselves it is instructive to describe problems which can 
take advantage of dual capacity of the computer, i.e. its analog/digital capacity 
in equipment type, organization, and concept. 

A problem in simulation can be well-suited to the hybrid computer in any 
of three ways : 

1. Because the physical system is best described in terms of both con- 
tinuous and discrete-valued variables. It is a hybrid system containing sampled 
data devices, logical control, or digital computers alongside equipments best des- 
cribed by transfer functions or differential equations relating continuous variables. 

2. Because the method of programming uses high speed continuous analog 
circuits - sampled and controlled from digital logic and storage equipment. The sys- 
tem might be a distributed system in which spatial variations of temperature, stress, 
potential, etc. is important to the simulation and thus a partial differential equa- 
tion must be solved by the iterative matching of boundary conditions. Within this 
class one must include cases where increased resolution requirements demand that an 
otherwise analog program must use digital eouipment for part of the simulation. 

3. Because the experimental procedure to be used with the simulation of 

a continuous analog-type physical system requires an information storage and evalua- 
tion procedure that is to be automatic. This represents a possibility for automatic 
design in which the design engineers evaluation of performance is replaced by a 
computer evaluation and parameter changes are effected automatically. 

Hybrid physical systems are becoming more common in aerospace design as 
more sampled-da ta , digital devices are incorporated. Typical devices are pulse radar 
systems, communications systems, on-off control systems, digital guidance computers, 
digital sensor. 

Hybrid programming, at present, relates to the use of analog circuitry 
(a few amplifiers, attenuators, multipliers, etc.) as a complex "arithmetic unit" 
with a digital program having iterative/sequential operation and information storage. 
Typically this technique is used to represent distributed flow and diffusion situa- 
tions as might be found in the design of ablative surfaces, rocket and jet engine 
nozzles, ion engines, nuclear reactor controls, etc. 

The second type of direct hybrid programming is well represented by the 
simulation of missile homing, rendezvous guidance, automatic landing systems, etc., 
where the transla tional motion and function generation is well simulated with digital 
equipment, and most other simulation recuirements are efficiently computed with ana- 
log equipment, providing the speed of simulation and output display features desired. 

Hybrid experimental procedures provide the analog simulator with the program 
and data adjustment capability of a stored program machine and promise interesting 
possibilities for the future. 
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THE EVOLUTION OF HYBRID COMPUTERS 



Hybrid computation came into being, as early as 1958, as a means 
of simulating the complete mission of certain aerospace vehicles. The 
task required computational speed exceeding the capabilities of the largest 
digital computers and a precision beyond that of the best analog computers. 
Combined systems of large general purpose analog and digital computers 
were created to solve this type of problem. The experience of these efforts 
and other experiments with analog and digital techniques in recent years 
have led to a growing hybrid computer technology. It is considered that these 
are the formative years in the evolution of a type of computer simulation in 
which the differences between analog and digital computing technique will 
dissolve through the development of general purpose scientific simulator 
system. 
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I. Some Historical Notes 

During the nineteen -fifties the capabilities of electronic computers 
expanded so fast as to stay well ahead of the needs of the average computer 
user. Such was the case in both the analog computer and scientific digital 
computer fields. One effect of this situation was the formation of two schools 
of experts with opposite views on the choice of the "best" general purpose 
scientific computer. Differences of training, experience, and semantics led 
to a serious barrier to communications between these two groups. 

The earliest attempts to combine the computation of analog and 
digital computers took place in about 1958 at the Convair Astronautics plant in 
San Diego and at the Space Technology Laboratories in Los Angeles. In both 
cases the job at hand was the complete mission simulation of the trajectory of 
a long range mis sile. The speed of the analog computer was a necessary 
element in the study to permit a "real time " simulation of the rapid motion of 
the vehicle and of control surfaces. However, the dynamic range required of 
the simulation was in excess of that of the best analog computers. That is, 
the ratio of the total range oJ^the trajectory to the required terminal phase reso- 
lution (a dynamic range of 10 to 10 ) was greater than 10 , the upper limit of 
analog computer dynamic range for small programs. Hence the digital computer 
was used to calculate those variables for which such dynamic range was neces- 
sary. The most important of these were the navigational coordinates -- the 
digital computer performed the open integration of velocities to determine the 
vehicle's position plus the dynamic pressure, a function that is very sensitive 
to altitude and velocity. 

It is fortunate that in such long range aerospace trajectory simulations 
the variables with wide dynamic range requiring precise calculation are not, 
at the same time, rapidly changing. Moreover the "high speed" variables do 
not require precise calculation. The early combined computer systems em- 
ployed the largest and fastest digital computers available at the time -- Univac 
1103A and IBM 704 -- together with 300 to 400 amplifiers of general purpose 
PACE analog equipment. In both cases even these fast digital computers were 
only just fast enough to perform the required repetitive calculations for the 
slowly changing variables of the simulation in real time. 

Since the installations of the first combined computers at least a dozen 
computer laboratories have employed general purpose analog and digital com- 
puters together to solve simulation problems, and a number of attempts have 
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been made to devise special purpose systems of analog and digital devices. 
Among the latter are the CADDA of the National Bureau of Standards and the 
"pulsed analog computer " of the MIT Electronic Systems Laboratory hybrid 
computers of a unique type are the combinations of a general purpose digital 
computer and a digital differential analyzer (dda), illustrated by the Bendix 
G15 with the DA-1 attachment and the Packard Bell PB2 50 with the Trice dda. 
The former system consists of a small, slow computer with an even smaller 
serial dda ($50, 000 and $100, 000 respectively). In contrast to this the Packard 
Bell system combines a small, medium speed computer $40,000) with a large 
serial -parallel dda ($500,000). Among the systems of general purpose com- 
puters, generally large analog computers have been combined with both large 
(IBM 7090) and small digital computers (PB 250, LGP 30). 

A brief analysis of the applications to which existing installations 
of combined systems have been applied leads to these generalizations. For 
the most part the analog computers in these systems have been employed in a 
normal manner to simulate the dynamic behavior of physical systems by solv- 
ing sets of non-linear, ordinary differential equations, while the digital com- 
puter has performed one or more of the following three functions: complex 
control logic, storage of arbitrary functions or sampled analog functions, and 
high precision arithmetic primarily for numerical integration. Examples of 
the applications are: 

A. Analog Computer Plus Digital Control Logic 

A system that in itself contains discrete control functions of 
continuous dynamic variables is appropriately simulated by a hybrid computer. 
The kinetics of a chemical process are simulated by continuous analog means 
while its digital control system is represented by a digital program. Similarly 
an aerospace vehicle with an on -board digital computer, control system, or 
autopilot is simulated by hybrid techniques. 

B. Analog Computer Plus Digital Memory 

A very common difficulty in the simulation of a chemical or 
nuclear reactor is providing an adequate representation of the transport of 
fluid in pipes from one point to another -- from reactor to heat exchanger. The 
simulation of this transport delay of a dynamic variable, such as the time varia- 
tion of the fluid temperature, is very nicely accomplished by the use of a digital 
computer for storage of the temperature function for a fixed, or variable, length 
of time. Digital computer memory has also been used effectively to store multi- 
variable arbitrary functions --an operation which is seriously limited in the 
analog computer. 
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C. Analog Computer Plus Digital Arithmetic 

This type of application is the "classic" one where the digital 
computer is used to perform precise, numerical integration of space vehicle 
velocities to keep track of the exact position of the vehicle over a very long 
range flight. 

It should be noted that a significant difference is apparent in 
the applications of computer systems employing a very small digital computer 
and those with large, very fast, computers. In general the small machines 
are limited to execution of control logic programs, one or two channels of 
transport delay simulation, or limited function generation programs. Since 
numerical integration and complex function generation by digital programs 
require considerable time for each calculation, for each discrete step in time, 
only the fastest digital machines can be used effectively for these tasks. 

II. The Elements of Hybrid Com puters 

A. Digital Computers 

Many conflicting factors influenced the choice of digital computers 
used in hybrid systems. Computer speed and economics have probably been 
dominant. Since there are so many computers on the market today that have 
sufficient speed and that span the complete range of prices, it is more instruc- 
tive to examine the features that are essential for hybrid computation. 

1. Speed. The speed of execution of arithmetic operations is most 
important, and this is a function of memory access time and multiplier speed. 
The access times of 2 to 5 microseconds are currently popular. This means the 
time for addition of two numbers is 4 to 10 microseconds. Multiplication and 
division take longer -- times of 15 to 40 microseconds are generally available 
and quite satisfactory. Overall program speed can be increased by the use of 
index registers -- three registers is desirable; more are useful. Special in- 
structions for subroutine entry, for executing commands out of sequence, and 
for testing and skipping can help increase computing speeds. 

2. Word Structure. The basic requirement is for a fixed point, binary 
word of at least 24 bits. Since round-off errors affect the last several bits a 
smaller word size would result in a dynamic range limitation of less the 10^ . A 
longer word may be useful in a few applications where fixed point scaling may be 
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difficult. Floating point computations may make things easier for the program- 
mer but should not be depended upon at the expense of computational speed. It 
may be noted that the equivalent of fixed point scaling is a necessary part of the 
analog program, and hence floating point operations may not prove as advantageous 
as for some all-digital programs. Decimal format and character oriented 
machines do not offer any advantages for hybrid computation, and usually they 
are slower than equivalent sized binary computers. 

3. Input/Output. High data rates in and out of core memory and any 
feature that minimizes loss of computing time for input/ output operations are 
highly desirable. In addition a fast, flexible means for communicating control 
signals to and from the analog section of the hybrid system is necessary. Three 
kinds of control signals are usually provided: interrupt lines, sense lines, and 
output control signals. It is by means of these controls that the sequential opera- 
tions of the digital machine are made compatible with the parallel simultaneous 
operations of the analog machine. Since communications must be made with 
many points in the analog computer a number of these control signals are needed. 
Interrupt signals from outside the computer stop the current sequence of cal- 
culations and force transfer to another sequence. Sense lines simply indicate to 
the digital program the current state of devices outside the computer which may 
be sensed by specific programmed instruction. Other programmed instructions 
will send control signals outside the computer on the output control lines. 

4. Memory. As noted above, the digital computer main memory should 
be a high-speed, magnetic core. Since most hybrid applications do not require 

a large memory for either program instructions or data, four, eight or twelve 
thousand words of core memory should suffice. Larger memories may be de- 
sired for special digital programs and larger hybrid problems when more ex- 
perience has been gained in this field, thus expandability of a memory to 16K words 
is a good feature. Newer computers are being introduced with small, very high 
speed, "scratch pad" memories. Such memories may have cycle times less than 
a microsecond and are used to store intermediate arithmetic results. This feature 
increases the overall computation speed of the computer. 

The normal manner of operating an analog computer involves a 
fair amount of non-computing time when the computer remains in the Hold or 
Reset- mode. These intervals may range from seconds to minutes while adjust- 
ments are made, pots are set, or recorders change, or while the programmer 
analyzes results. It is not possible for the analog computer to operate on other 
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programs at these times, however, with a hybrid system, where such waiting 
periods are likely to occur also, it is reasonable to consider having the digital 
computer work on a different program during the intervals, whatever their 
length. Appropriate "interrupt" and "memory lock-out" features are pos- 
sible to permit time sharing of the digital machine without affecting the hybrid 
program and without the danger of one program interfering with the other. 
The secondary program (a strictly digital problem) is simply stored in a 
"protected" part of the core memory and utilizes all the bits of time not re- 
quired by the hybrid program. 

5. Peripheral Equipment. In many digital computer installations 
the investment in peripheral equipment rivals that in the central computer. 
Current hybrid computer applications require only a minimum of digital peri- 
pheral equipment. The graphic output equipment associated with the analog 
computer is sufficient for computational results. Punched paper tape reader 
and punch and typewriter may be all that is required for programming. Larger 
systems in the future will employ punched cards and magnetic tape, primarily 
for rapid change-over of problem and automatic check-out. Large off-line 
data storage does not appear necessary for most applications. 

In summary the digital computer must be characterized as a je- 
quential machine (Fig. 1). For effective use within a hybrid system the 
machine (a) must have sufficiently high internal speed for it to appear as though 
a number of calculations were taking place simultaneously; (b) must be or- 
ganized for maximum speed in executing mathematical calculations and (c) 
must have efficient means for input and output of data during calculation. 

B. Analog Computers 

In contrast to the above the analog computer is a para l le l machine 
with many computing components and I/O devices operating in concert. There 
are few, if any, features of the modern analog computer that are not appropriate 
to a hybrid system. However only the largest analog machines have be em used 
for general purpose hybrid simulation. The common measure of a large com- 
puter is that it has 100 to 200 operational amplifiers. Since two or more com- 
puters may be "slaved" together, larger systems are possible when required. 

Analog computer features that are important for hybrid systems 
can be simply listed as: 
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FIG. I ORGANIZATION OF DIGITAL COMPUTER 
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Integrators with multiple time scales 

Amplifiers for tracking and storing voltages 

— Very fast control of the modes of individual amplifiers 

Automatic, remote control of the setting of potentiometers 

Fast, accurate multipliers and trigonometric resolvers 

High speed comparators with logic signal outputs 

Electronic switches (logic signal gating of analog signals) 

In the early days logic equations or switching functions were programmed with 
relays and stepping switches, which were connected to the patch board by 
various means. Present day technology employs electronic switching of in- 
tegrator modes and voltage signals at high speeds, and the delays inherent in 
relay devices can no longer be tolerated for logic operations. The logic build- 
ing blocks common to the digital computer designer (flip-flops, gates, inver- 
ters, monostable multivibrators, shift registers, and counters) are ideally 
suited to these operations. Thus with electronic switches replacing relay con- 
tacts, logic modules have become an integral part of all new, large, analog 
systems. These modules are programmed like the other analog components 
by interconnections at a patch panel. Many signals occur simultaneously but 
they are logic signals -- two values, Zero and One, that change as functions of 
time. Input signals to logic programs come from comparators, manual switches, 
and external control signals. Logic program outputs go to integrator mode con- 
trols, storage amplifier controls, electronic switches (DA switches) to gate 
analog signals. As will be shown later, it is essential for a hybrid system to 
have a very significant complement of digital logic components. A few gates 
and flip-flops are not sufficient. The requirements for use of logic components 
in an analog computer for hybrid operation are so great that the EAI HYDAC Digital 
Operations System is an entire computer console with its own patching system 
used entirely for the programming of digital components for parallel computation. 
This console is really a complete logic computer . It is used together with a 
conventional analog computer to form what is truly an all parallel hybrid computer. 

In summary, the modern analog computer must be characterized as a 
parallel machine (Fig. 2). It is not solely for continuous variables. It is a 
parallel assemblage of building blocks: integrators, multipliers, etc. , for con- 
tinuous variables. It is organized for convenient representation of an "analogous" 
physical system by means of a computer model constructed of these building 
blocks. 
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FIG. 2 ORGANIZATION OF ANALOG COMPUTER 
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C. Conversion. Devices 

In providing data communications between a sequential computer 
and a parallel computer, three kinds of devices are commonly used: the multi- 
plexer, the analog -to-digital converter (ADC), and the digital-to^analog converter 
(DAC). In addition, all early systems have employed a timing and control unit 
which performs a relatively fixed set of operations, with manual switches 
to select options such as sampling frequency, and number of channels. Such 
"linkage systems" thus consisted of a timer unit plus a group of linkage 
building blocks prewired to perform a specific task. With the integration of 
digital logic components into the parallel computer, however, greater program- 
ming flexibility is possible by use of these logic units for timing control of the 
data conversions. Furthermore the converters and multiplexer can act very 
naturally as additional building blocks in the parallel computer. Thus it is likely 
that future hybrid systems will simply incorporate the "linkage system" within 
the parallel computer. 

Usually several or many analog signals in a hybrid program will 
be sampled, converted and transmitted periodically to the sequential digital 
program. The numbers, of the several sequences of numbers to be entered 
into the core memory, can be accepted only one at a time. Since this is so, the 
conversions from voltage to number form can be performed one at a time -- 
first from one analog variable and then from another. The multiplexer is used 
to select one from many analog signals, to step through a sequence of signals, 
and thus to furnish voltage input signals to the ADC. 

The output of commonly used ADC's is a binary number of 10 to 
14 bits. A 13 bit binary output probably is the best compromise; for it represents 
a resolution of one part in eight thousand, and resolution of analog voltage sig- 
nals is at best one part in ten thousand. Conversion times range from 50 to 300 
microseconds. A typical time of 100 microseconds would allow the converter to 
be shared by 16 analog signals each with a frequency spectrum extending to 20 
or 30 cycles per second. 

DAC units should lave the same binary word size as the ADC, 
except for special low accuracy uses. Conversion times are not determined by 
the converter but rather by the bandwidth of the analog amplifier following the con- 
verter. Data from the sequential computer appear only one word at a time, and 
some means for retaining the latest value, of each sequence of numbers, for each 
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output function is needed. The sequence of numbers coming from the computer 
may first be converted to voltage values by a single DAC, and then distributed 
to storage amplifiers for each channel. It is more customary, however, to 
hold each of the latest words for each channel in a digital register which is 
an integral part of the DAC assigned to each channel. (Fig. 3) 

D. Special Forms 

As a passing thought it may be noted that while the primary em- 
phasis here is being placed on the distinction between parallel and sequential 
operation, the term "hybrid", historically, has been used to imply the combi- 
nation of continuous and discrete calculations, and that therefore consideration 
might be given to two special kinds of hybrid computers: 

1. The Parallel Hybrid Computer, which is a proper term for the 
EAI HYDAC 2000 machine. This system combines an analog computer with a 
general purpose system of programmable logic building blocks, multiplexer, 
ADC, DAC's digital memory units for storage of sampled analog functions, and 
several digital numerical adders and subtractors. The application of this system 
encompasses an intermediate range of hybrid problems, such as: 

a. Transport delay simulation 

b. Single and multi -variable function generation 

c. Logic control systems 

d. Automation of the analog computer for parameter searches 
and optimization studies 

e. Simulation of numerical and sampled data control systems 

2. The Sequential Hybrid Computer, which is exemplified by the ex- 
perimental "pulsed analog computer" techniques developed at M.I. T. for use in 
an aircraft flight trainer. This system employs a sequential digital computer 
which controls a small number of analog functional components -- one multiplier, 
one reciprocal generator, one integrator, and several adding units. These units 
are interconnected and receive inputs by digital program control. They form, 

in effect, "analog subroutines" for the sequential computer. 

Ill . The Se quenti al Par allel H ybrid Compu ter 

A. Description of System 

The term "hybrid" is most appropriately used to indicate the 
combined use of sequential and parallel computing techniques, first because the 
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future growth in hybrid simulation will be predominantly in this direction. 
Second, and more important, is that from the standpoint of the programmer 
who must bring the two types of computers together to find a useful solution 
to a problem, the only really significant disparity lies between sequential and 
parallel operations. The difference between continuous voltage and discrete 
number is simply one of format. It would make little difference if the analog 
signals were frequency modulated, pulse code modulated, hydraulic, or 
pneumatic -- appropriate format converters could be found. The feature of 
hybrid computation that is of importance is: in part of the machine many 
operations are taking place simultaneously, and many, time-varying problem 
variables exist in parallel; while elsewhere a number of operations take place, 
one at a time in a repetitive manner, so as to effect the generation of several 
problem variables, as if they occurred simultaneously. Furthermore, the 
parallel computation is tied to a real time base: the very passage of time 
itself accounts for the changing of the basic independent computer variable. 
The sequential program is asynchronous -- not controlled by a clock. Opera- 
tions are executed in sequence at whatever rate is possible, and for any refer- 
ence to be made to the actual elapsed time external communication is necessary. 
This basic incompatibility requires that the interface between the two types of 
operation embody more than the simple format conversions performed by the 
ADC and DAC's. It is necessary for data and control information in the paral- 
lel machine to be available to the sequential machine and conversely that the 
latter be able to send data and control signals to many points in the former. 
Coincidence or simultaneity of events communicated to and controlled from the 
sequential program are particularly difficult to handle. The logic and data 
control of the interface equipment must resolve these differences in timing and 
operation. What might be termed an "impedance matching" device is needed 
between the parallel and sequential program in order to make most efficient use 
of both machines. The exact manner in which this is done will vary from problem 
to problem. (Fig. 4) 

B. A Simple Example 

An example will illustrate some basic considerations in defining a 
general purpose hybrid system. First the operation of the simplest of linkage 
systems: an ADC and multiplexer, a number of DAC , and an "interrupt clock". 
The flexibility of the stored program digital computer is relied on for control of 
these units. Assume 10 analog signals are to be converted one at a time. These 
words are placed in memory (average program time: 40usec. per word) and then 
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about 7. 5 or 8. 0 milliseconds of sequential, digital calculations takes place, 
followed by output from memory of ten words (20 sec. per word) to ten DAC's. 
The entire cycle requires 7500 +10 x 40 (input) + 10 x 20 (output) - 8100 micro- 
seconds of digital program time. If it is assumed the conversion of the data 
( A to D ) requires 150 microseconds per word then 1. 5 additional milliseconds, 
or 9 ms are needed if everything proceeds sequentially. Assume further that 
because of the frequencies of the analog signals it is necessary to sample at 
least some, and therefore all, of the channels at 100 samples per second. A 
"real time interrupt clock" is set at 100 cycles per second. This timer unit 
is an adjustable oscillator that sends an interrupt pulse to the digital computer. 
The latter then activates the ADC, waits 150 m seconds for a completion, or 
"Ready" signal, steps the multiplexer to the next channel, stores the converted 
word in memory, and then repeats this cycle ten times. With the tenth step 
the multiplexer resets to the first channel. The program then proceeds with 
the 7. 5 milliseconds of calculation, outputs ten words one at a time to ten DAC's 
and then waits for the next interrupt pulse. (See Fig. 5, 6, and 7. ) Manual controls 
are provided for selecting the interrupt clock frequency and the number of chan- 
nels in the multiplexer stepping cycle. 

C. Problem Areas 

This is certainly a simple system and it appears to satisfy the 
basic requirements for communications. Some of the shortcomings of the system 
are apparent: Sampling and outputting of each channel do not take place simul- 
taneously; 15% of the sequential program time is "waiting" time, and 3 to 5% 
is used to select and control devices external to the sequential computer. In 
other problems these percentages may be higher. The other weaknesses in this 
system lie in the fact that it was not designed to be a general purpose system. 
For example, the placing of logic, timing, control and storage functions in the 
digital computer severely limits its arithmetic processing capability. This ap- 
proach is restricted in application to a class of problems for which the periodic 
"input/ calculate/ output " cycle is useful. 

D. System Improvements or Variations 

By programming the parallel digital components of the parallel 
computer to perform timing and control functions for the system the following 
changes to the above system are suggested: 

1. Simultaneous sa mpling . If the sequential program operates on two 
or more of the input numbers together to calculate an output, then errors may 
occur since the input numbers were sampled at different times and correspond to 
different values of the independent variable. A similar effect may occur at the 
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output since the numbers in a group of ten appear at the ten DAC's at different 
times. It is certainly possible by numerical means to compensate for the er- 
rors, at the penalty of additional program time. The common solution is to 
add memory to eac.h of the ten input and ten output channels. Ten Track/Store 
amplifiers are added before the multiplexer and a control signal causes them 
all to sample, by storing the voltages, simultaneously. At the output, 13 bit 
registers are added in front of the DAC's. When all ten registers have been 
loaded, a transfer pulse causes all DAC values to be updated at the same time. 

2. External timing of ADC and multiplexer . Sequential program time 
can be saved by permitting the control of the ADC, Multiplexer, TrackStore 
cycle to be controlled externally. (Fig. 8). A simple clock, counter, flip-flop, 
and group of logic gates will permit the input conversion cycle to run at its 

own rate -- interrupting the sequential program only at the completion of a 
conversion. Thus the conversion time can overlap the calculation time, elimina- 
ting the waiting time. Upon interrupt only 10 to 20 seconds may be required per 
sample; many control steps are eliminated. Similarly on output, the addressing 
and selection of output channels can be done by simple circuits rather than using 
sequential program time. 

3. Real t'ime clock to est ablish sa mpling frequency . If the se- 
quential calculation involves numerical integration over a long term, the ac- 
curacy of the sampling interval is just as important as the round -off and trunca- 
tion error in the numerical calculation. Although numerical means may be resor- 
ted to for very accurate integration, in a hybrid program the calculations still 
need to be referred to a real time base. This is done by using an accurately cali- 
brated source for setting the sampling interval, or frequency. A good, high- 
frequency, crystal stabilized oscillator is an important part of the parallel digital 
subsystem. Sampling frequencies lower than the oscillator frequency are 
selected by use of preset counters. 

4. U se of sense lines to reduce number of conversions . In the simple 
example problem only whole number data are transmitted to the sequential pro- 
gram . Thus, if the relative magnitude (greater or less than) of two analog 
signals is needed in the di gital calculation, the two numbers must be converted, 
stored and then compared. This can be accomplished more simply by use of an 
analog comparator the output of which is sent to the digital computer by a sense 
line -- saving time and equipment. The state of any parallel logic component may 
be monitored conveniently by sense lines. These are tested in one memory cycle 
(2 - f}i(sec. ). If many such communications are needed the savings will be sig- 
nificant. Sense lines should also be added to allow the sequential program to 
determine the modes of the analog computer, the relative sizes and signs of error 
quantities, and the states of recording devices. 
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5. D etection of random ev ents. With fixed, periodic sampling 

the sequential program cannot tell exactly when events take place in the paral- 
lel machine. With comparators and parallel logic, complex functions of analog 
variables can be monitored. For example, it might be required to determine 
when the overshoot in x^ exceeds x^ after the third cycle; but only when x^ 
is negative and x is less than x^. After determination the sequential program 
can be interrupted to perform specific conversions and calculations -- 
asynchronously with respect to the primary conversion cycle (Fig. 9). In 
this manner the parallel logic avoids the delays in the sequential program 
and uses the latter only when required. The parallel logic program analyzes 
the data, interrupts the sequential program, sets up the proper channels for 
conversion in and out of the digital computer. 

6. Multiple samplin g fre quencies. In the example problem all 
channels are sampled at a frequency determined by the highest frequency present 
in any one channel. It may often be the case that there are two or more groups 
of variables with different ranges of variable frequencies. It may then be ap- 
propriate to sample each group at different frequencies. Another approach 
using different sampling rates is to use several eight channel multiplexers in 
cascade so that the output of two of them feed two channels of a third which 
feeds the ADC. On each cycle of the third unit the first twc^are stepped, yield- 
ing different variables for those two channels on each cycle. Alternatively, 
each time the third steps to the two special inputs the corresponding multiplexer 
makes a full cycle. Timing control of these operations is performed by paral- 
lel logic components. (Fig. 10). ^ ^ 

7. Asynchronous sampling. A completely asynchronous conversion 
system has been designed by one computer laboratory, in which the sequential 
program is interrupted only by comparators. Twenty analog problem variables 
are compared to reference values that are adjusted by the digital computer when 
necessary. Each comparator calls for conversion of some group of variables 
(the same variables may be called for by different comparators). When two or 
more comparator signals occur simultaneously or during a conversion operation, 
two levels of priorities are set up by logic elements to determine what interrupts 
are to be made. While the system appears complex, it is accomplished in a 
simple fashion in the parallel computer and makes good use of the sequential 
computer time. 

A longer list of useful variations in the control and timing of se- 
quential/parallel communications can be compiled. It can be said, however, 
that the omission of the above capabilities in the hybrid computer severely limits 
the flexibility and high speed capability of the computing system. For the most 
part, these features should be explained in terms of the particular problem ap- 
plications. 




COMP 



> 



FIG. 9 DETECTION OF RANDOM EVENTS 



Hybrid Seminar 
Page 22 



E. Operating Times for Typical Mathematical Functions 

The repeated emphasis upon the efficient utilization of the 
sequential program time, high arithmetic speeds, and programming tricks 
to gain speed can be seen to be warranted when one examines the sequential 
operating times for several typical mathematical functions, which, on the 
analog computer, would be executed continuously and in parallel. 



a. The sum: 



a + b+ c + d 



40 microsec. 



b. The expression: 



ax+ by + cz 



160 microsec, 



c. Sin«t or cos art 

2 2 

d, Square root of x + y 



215 microsec 



432 microsec 



e. Generate z - f (x,y), where two 

dimensional interpolation is required 
between functional values evenly 
placed in x and y: 



0. 5 to 1. 5 millisec. 



f. Rotate a vector through three co- 
ordinate angles: 



2 to 6 millisec, 



Perform one integration of a single 
derivative for a single time step: 



0. 1 to 1. 3 millisec. 



Allowing another millisecond or two for control and input /output instructions 
one can estimate the real time speed performance of this program. The speed 
is best expressed in terms of the useful upper frequency (at full scale) in a 
problem variable that can be represented by the computer . Although the exampli 
equations have no real meaning the frequency limit for such a program is about 
11/2 cps. This does not seem like very fast performance for such a few simple 
equations. On the other hand it is fast compared to frequencies of some of the 
variables in an aerospace simulation program for which digital precision is re- 
quired. In any event it should be clear for best utilization of the sequential 
computer care should be taken to reduce operating time whenever possible. 



IV. HYP AC 2400 System Description 



Experience derived from the successes and failures of the early 
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hybrid systems as well as considerations such as those above led to the first 
general-purpose hybrid computer, the HYDAC 2400. This computer satisfies 
all prerequisites of a compatible system of parallel and sequential computing 
components as discussed earlier. 

In order to facilitate the understanding of the programming aspects 
of an aerospace problem on the HYDAC 2400 discussed in the next chapter, a 
brief description of the system is given here. 

The HYDAC 2400 computer system is a successful integration of 
general -purpose analog and digital computers. It consists of three computers: 
the 231R-V analog computer, the DOS 350, and the 375 ( 3C DDP-24) digital 
computer. (For system block diagram, see Fig. 11) In the following, a short 
outline is given of each section of the HYDAC 2400, and the communication be- 
tween sections. 

A. 231R-V Analog Computer System 

The 231R-V analog computer consists of standard computing com- 
ponents under control of a versatile, highly sophisticated control system. The 
heart of this control system is the analog memory and logic system (MLG). The 
MLG system adds a capability for programmed multiple time base and multi- 
speed repetitive and iterative computer operation. The system utilizes electronic 
mode control and high-speed analog memory in order to maintain the accuracy 
necessary for "real-time" operation. Various forms of analog memory (micro- 
store, track and store amplifiers) are essential for simple and efficient data 
transfer into the analog computer. The separate patchboard allows the controls 
to be pre-patched, or control commands from external digital computing equip- 
ment may be introduced to vary the program during problem solution. These 
control signals may affect the operation of the analog computer or some of its 
components many different ways: 

(1) change mode of operation 

(2) change time constant of an integrator 

(3) control analog memory 

(4) set or monitor a potentiometer , or other components 

(5) make a logical decision on an AND gate, etc. 

The high accuracy solid state multipliers, resolvers and function generators are 
also important contributors to the overall system performance. 
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B. 375 Digital Computing System 

The 37 5 forms the arithmetic and primary memory 
capability of the HYDAC 2400 system. It is normally controlled by the DOS 
350. 

The 375 is a high-speed, all solid state general purpose 
digital computer with 24-bit word length and sign -magnitude binary number 
representation. The size of the random access memory can range from 4,096 
words to 32,768 words, with a memory cycle time of 5 microseconds. The 
lO^sec. addition and 31^.sec. multiplication times are typical of the speed of the 
arithmetic unit. 

Standard input/output capabilities include four data input/ 
output channels: A buffered character input channel, a buffered character 
output channel, a parallel 24-bit input channel and a parallel 24 -bit output 
channel. The automatic interrupt may be used in connection with any channel 
desired. In addition,l6 simple sense lines and 16 output control lines are 
provided for efficient control signal communication as will be described later. 

C. Ihe DOS 350 Digital Operations System 

The DOS 350 is the control and communication center for 
the HYDAC 2400. Its modular, functional structure of digital computing and 
conversion components satisfies the stringent requirements established by the 
basic differences of the parallel analog and the sequential digital computers. 

The logic building blocks, counters and arithmetic elements 
provide an extremely high speed parallel counting, decision making, control and 
arithmetic capability that cannot be satisfactorily provided elsewhere in the 
system. The same is true for the optional storage units that may be included 
in the DOS 350. 

The communication system provides for the flow of infor- 
mation throughout the system. The information handled can generally be classi- 
fied as either whole word data or as bits or signals that must be communicated 
between the sections. The conversion of whole word data is performed by the 
analog -to -digital and digital -to -analog converters. The appropriate channels 
for conversion are selected and sampled by the multiplex equipment. 

In addition to the transfer of whole word data, provision must 
be made for the communication of logic and control signals. The sense lines 
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enable the DOS 350 to exercise control over the digital program. The 37 5, 
through the use of a specific test instruction (SKS) is able to determine their 
status and make consequent decisions. When the DOS must have immediate 
access to the 375, an interrupt line has to be activated. When this happens, 
the 375 is forced to store its present computation and proceed to a program 
designed to handle the interrupt situation. The system has up to eight inter- 
rupt lines each of which has a specified priority in relation to the others. 
Similarly, the 375 also has the ability to transmit signals to the digital section. 
Output control pulses are generated by the 375 and terminated on DOS 350 
patch panel. 

The control signal exchange between the 231 R-V and DOS 350 
computers is done by means of special interface components: these are the 
analog comparator and the digital analog (DA) switch. The analog comparator 
permits the DOS 350 to be made aware of the occurrence of certain events 
within the analog section. This circuit generates a logical ONE on the DOS 350 
patchboard when the algebraic sum of two selected input voltages exceeds zero. 
The DA switch is basically an analog gate; under DOS control, the path of an 
analog signal may be open or /closed. Thus the logical program is able to 
accomplish numerous functions, such as changing parameter s , providing 
automatic rescalings, etc. 

While the above description of the HYDAC 2400 is extremely 
brief, it is hoped that it would aid the reader in understanding subsequent dis- 
cussions. More detailed writings on the systems are available from numerous 
sources. 
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HYBRID SIMULATION OF A RE-ENTRY PROBLEM 



During the past ten years a great deal of research and development 
work has been conducted on various types of re-entry vehicles. Numerous tech- 
niques for guiding and controlling re-entry vehicles have been proposed. The 
purpose of this simulation is to evaluate the capability of a new flight control 
system for re-entry vehicles which by its nature is simple, reliable and 
inherently insures a safe re-entry. The ensuing discussion explains the 
general re-entry problem, the unique control system being studied, and includes 
a detailed discussion of the simulation equipment required and its programming. 

I. General Re-entry Considerations 

To obtain a better insight into the need for a simple yet reliable 
guidance and control system for lifting re-entry vehicles, a brief discussion 
of the general re-entry problem follows. 

The re-entry vehicle utilized in the study is an unpowejred lifting 
vehicle with wings, which, unlike its ballistic counterpart, is highly maneuver able 
and can be landed on conventional runways. (See Fig. 1) 

The lift and drag coefficients (C L and ) of a typical high L/D 
lifting re-entry vehicle are shown in Fig. 2 as a function of angle of attack. 
A plot of lift to drag ratio (L/D or C /C-q) * s also shown since it is of importance 
in determining the range of the re-entry vehicle. During the high velocity portion 
of re-entry flight the vehicle operates on the high C (large*?) side of (L/D) max. 
(15°< < 65 ). These larger lift coefficients yield a re -entry trajectory with 
lower dynamic pressure, acceleration, and temperatures. Fig. 3 shows a typical 
uncontrolled re-entry trajectory with an angle of attack, of 35 . Fig. 4 defines 
the coordinate system and symbols. During most of the vehicle flight r e *V*^ tr 0 
and)fe?0. Thus, the following approximate equations can be used: 

,4k tKA-th- vet, f 5 , 
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For a lifting re-entry vehicle the lift is large enough to allow the vehicle to 
fly along an equilibrium glide path where hs^O. Along the equilibrium glide 
path where ha* *0 , 



where q is the dynamic pressure PV ) and S is the wing area. Thus for 
a given there is a unique altitude -velocity profile. Fig. 5 shows the altitude - 
velocity profile for two different equilibrium glide lines and the boundaries of 
the re-entry corridor. The lower boundary of the corridor is determined by 
the temperature and load limits of the vehicle and the upper boundary by the 
recovery ceiling. The recovery ceiling is defined as that altitude (with 
altitude rate, h«0) from which the vehicle can recover without exceeding 
the temperature and load limits of the vehicle. 

One of the primary re-entry problems is to control the vehicle 
so that temperature and load limits are not exceeded and a smooth equilibrium 
glide is established. The heavy trajectory shown on Fig. 5 is a typical uncontrolled 
re-entry with its familiar skipping oscillations which cause the vehicle to approach 
dangerously close to the heat limits. The temperature rate control system being 
simulated, was developed to eliminate these skipping oscillations and reduce 
the peak temperatures during re-entry. 

Another re-entry problem is to carefully manage the energy of 
the vehicle as it re-enters so that the desired terminal point is reached. The 
range capability of the re-entry vehicle can be determined quite readily in 
the re-entry portion of flight since the vehicle is near equilibrium glide where 
hftjO. If both sides of Equation 2 are divided by L/M or Q - the following 

is obtained: 



V _ _ D ( 4) 
dV 

If V "Jj^ is substituted for V: 

dR=(u/D)^l (5) 
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If this is integrated from an initial velocity to zero, the range is obtained: 



x n>9 7 

Thus, daring the re-entry portion of flight the range is a function of velocity 
and L/D. Range control is thus accomplished by varying L/D, i.e. , by varying 
angle of attack. From Equation 6 it should also be noted that the range is 
very sen sitive t o initial velocity when the velocity is nearly equal to orbital 
velocity ^ t© g . For the re-entry shown in Fig. 3 the sensitivity of range 
to initial velocity error is approximately 80 NM/fps. 

Lateral maneuverability is obtained by banking the re-entry 
vehicle so that the aerodynamic lift vector is rotated providing a lateral 
acceleration. Fig. 6 shows an energy management footprint for a typical 
re-entry flight. The lines of constant <K and show what attitude must be 
maintained to reach a particular landing site. The dashed lines show temperature 
limits. This large maneuverability of a lifting re-entry vehicle requires a 
reliable guidance system which will perform accurately over the long re-entry 
and minimize errors at the desired terminal point. The temperature rate 
control system is being simulated to demonstrate its compatability with different 
types of navigation and guidance systems. As will be pointed out in the next 
sections the TRFCS acts as a filter to the guidance signals to insure the safety 
of the vehicle at all times. 

II. The Statements of the Problem 

A. Problem Background 

The temperature rate flight control system (TRFCS), developed 
by the AC Spark Plug Division of the General Motors Corporation, is based 
upon the use of temperature sensors instead of conventional inertial instruments 
to provide both short-period stabilization and long-term guidance during the 
re-entry flight. Details of this flight control system are given in the accompanying 
AC Spark Plug Report. (See Reference #1) . The mathematical formulation for 
the simulation of the re-entry problem was furnished by AC Spark Plug. 
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The new control system introduces several significant advantages: 

(1) Overall vehicle safety during re-entry. This system repre- 
sents an unorthodox approach to the design of an overall 
system of re-entry vehicles. In the standard approach the 
tendency has been towards a complex integrated system, 
while in the TRFCS, successful effort has been made to 
separate safety of the vehicle from the task of accurate 
navigation. Because of the inherent nature of temperature 
rate feedback and certain selected limits on the control 
authority, the control system minimizes skin temperature 
peaks. The maximum "g's" amd dynamic pressure are 
independent of initial conditions and maneuvers performed. 
This safety aspect of the TRFCS performance is entirely 
independent of the guidance commands and in fact, the TRFCS 
serves essentially as a filter for them. 

(2) Simple, reliable hardware. This separation of control and 
guidance also results in more reliable hardware, since the 
failure of the necessarily complex guidance system cannot 
cause the complete destruction of the vehicle. Furthermore, 
simple thermocouple temperature sensors replace the con- 
ventional gyros and accelerometers. These sensors are 
used to control the flight path as well as the short -period 
oscillations in pitch and yaw. The only additional sensor 
required, besides the temperature sensors, is a vertical 
reference gyro, which, for the safety aspects of the re-entry, 
can be quite inaccurate. 

(3) Both manual and automatic modes. In case of automatic 
guidance system failure, the TRFCS can be controlled 
manually. The manual flight program to be followed by the 
pilot, is very simple and the resulting temperature peaks, 

dynamic pressure, and "g n loads compare favorably 

with those obtained in the fully automatic mode. 

During past years, extensive simulation studies were conducted by 
AC Spark Plug to obtain familiarity with the control system. A rather conventional 
simulation program was pursued: First analog simulations were performed to 
gain qualitative knowledge of the system, and to determine the practibility of this 
approach. Next, digital techniques were used to evaluate the accuracy of the 
guidance through TRFCS. 

In the analog simulations, the system characteristics were split 
and analyzed in two independent studies. In the first, a three degree of freedom 
simulation of the mass center of the vehicle was combined with equations des- 
cribing the short-period pitch dynamics of the vehicle. Pitch axis controls and 
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trajectory controls in three dimensions included an approximate, simple lag repre- 
sentation of the lateral response of the vehicle. In the second type of simula- 
tion, the effects of lateral dynamics were obtained by simulating the dynamics 
of the vehicle in detail by a standard set of lateral stability equations with 
variable coefficients. The coefficients were varied with dynamic pressure, 
velocityjand stagnation temperatures of the vehicle skin, all obtained from 
function generators. The data for setting up the function generators came 
from the first type of simulation. In turn, the results from the second type 
of simulation were used to determine the lumped lateral response for the first 
type of simulation. Thus, a basis for an iterative procedure was established. 

The reason for separating the simulation of the pitch dynamics 
and trajectory control from simulation of the lateral dynamics was the un- 
availability of necessary simulation equipment. The investigation of aspects 
of the system such as coupling between pitch and roll, could not be made with 
the above "split" simulation approach and awaits a complete six-degree of freedom 
simulation. 

The next logical step was therefore to study the system's 
characteristics in a combined, six-degree of freedom simulation. But here 
arises the question of what computer or computers should be used. Past 
experience has shown that the conventional first -analog then digital approach 
is definitely not the best approach. Some of the conclusions gathered during 
the simulations are as follows: 

(1) Repeatability of analog simulation was only marginal 
(50 miles in range). 

(2) Slowness of digital simulation. Even for narrow ranges, 
determined by previous analog simulation, digital simulation 
was too time consuming, and therefore too expensive to 
optimize parameters. Reduction of the digital data also 
proved to be a problem. 

B. Problem Objective 

The main objective of the problem is to evaluate a TRFCS 
controlled re-entry both in automatic and manual modes of operation. To fully 
explore the ability of the TRFCS to control the short period attitude of the 
vehicle throughout the re-entry, a complete six degree of freedom simulation 
is required. To evaluate the ability of the TRFCS to guide the vehicle during 
re-entry to the desired terminal pointy an accurate and repeatable simulation 
is required. Economy of analysis should be considered, especially in the 
automatic guidance studies where faster than real time simulation can be 
employed. 
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C. Computational Requirements 

In order to attain the above problem objective, the following 
set of rigid computational requirements must be met: 

(1) High accuracy in trajectory calculations for the evaluation 
of the guidance capability of the TRFCS. 

(2) Very fast computing capability to faithfully simulate the 
high frequency parameters for the short -period dynamics 
of the vehicle. 

(3) Real time and faster than real time simulation for control 
system evaluation. (For economical evaluation of the 
control system in automatic mode, the time scale should 
be as high as possible,). 

On the basis of experience gained during past simulations, it was 
concluded that the simultaneous need for high accuracy and very fast computation 
can only be satisfied by a hybrid digital -analog computer. 

Such a computer would allow the programmer to choose either 
analog or digital solution for different portions of the problem, trading it with 
fast processing for high resolution etc. 

A hybrid system makes high demands on both the digital and the 
analog elements as well as on the control section. The tasks to be performed by the 
control section, however, cannot be over emphasized. Some of these tasks are 
the basic timing, control operations, logic decisions and conversions which 
preferably are parallel with other computations. Only in this way can the pro- 
grammer truly utilize all elements of the hybrid to the fullest. If the digital 
section for "useful" computation is severely limited. To achieve these complex 
tasks, the "control center" should consist of programmable logic elements, such 
as flip-flops, counters, shift registers, parallel memory and converters. The 
answer to all these requirements was EAI's HYDAC 2400 system. The selection 
of the HYDAC system enabled the programmers to use analog and digital compu- 
tational equipment judiciously and to produce the best mechanization of the problem 
in question. 
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II. Problem Mechanization 

All advantages of the hybrid computer are in vain, unless very 
careful consideration is given to the programming of the physical system under 
study. This phase of the simulation transforms the essentially general purpose 
computer into a simulator of the specific physical system. 

A. Allocation of Tasks 

The first step towards a successful hybrid program is the allocation 
of tasks on the computer. The underlying philosophy is to subdivide the physical 
system into sections, and assign these to various parts of the computer, where 
their speed and accuracy needs are best satisfied. 

As shown in Fig. 7, the physical system consists of four sections, 
three of them (the vehicle dynamics, the TRFCS, and the temperature sensor 
simulation) constitute the attitude control loop, while the vehicle dynamics, together 
with the guidance system and TRFCS, form the long period loop. 

The assignment of these sections to various elements of the 
HYDAC 2400 is shown in Fig. 8. The attitude control loop, consisting of the 
vehicle rotational dynamics, the TRFCS and the short period sensor equations, 
are programmed on the analog section. In addition, the displays and cockpit 
simulator are tied into the analog since continuous analog signals are required. 
The translational equations of motion, long period heat sensor equations and 
guidance equations are programmed on the 375 because of the stringent accuracy 
requirement. The DOS 350 provides the master timing, data conversion, function 
generation, and reaction jet control logic. 

The DOS 350 timing and control is essential because of the 
operational differences between the analog and digital sections. The analog is 
a parallel continuous computer with computing time independent of problem 
size. The digital is a serial, discrete inter -val computer with computing time 
directly dependent on the size of the problem. The DOS, through its timing and 
controls^ synchronizes the calculations on each computer and controls the flow of 
information between sections. Function generation and the reaction control jet 
logic are ideally suited to the DOS since these operations can be performed rapidly 
in parallel with the 375 -general purpose digital computer, so that the digital 
computation time is minimized. 
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B. DOS 350 Program 

Fig. 9 shows a detailed block diagram of the DOS. The DOS 
performs the five following functions: 1. ) Timing 2. ) Mode Control, 3. ) Data 
Transfer, 4. ) Function Generation, and 5. ) Reaction jet control logic. These 
functions are described in detail in the following sections. 

1. Timing 

Timing is required on a hybrid computer for the following reasons: 

(1) To make the mathematical time step used in the integration in 
the digital program correspond to the physical time scale used 
on the analog computer; i. e. , the digital and analog sections 
must run in synchronism. This synchronism is accomplished 
by sending a periodic master time pulse T^ which initiates the 
calculations for each time step in the digital computer. 

(2) To time information transfers between the analog section and 
the digital section. Not all transfers are at the same rate, 
since the serial memories of the DOS are used for function 
generation of aerodynamic moment and force coefficients. 
These variables must be transferred at a high rate, since the 
functions are used in the short period rotational dynamics of 
the system. On the other hand, those variables relating to 
long term trajectory variations are transferred to the digital 
section at a lower sampling rate. Timing is clearly needed 
to control these two different transfer rates. 

Timing is controlled from the DOS with a Master Timer which 
controls the sequence of events occurring per cycle of operation. The master ^$0 
timer is a BCD counter driven by the SM8 signals occurring at the rate of ^^pps 
(approximately one pulse every 2 ms. ) Fig. 10 shows a block diagram of the 
master timer. 

As shown schematically in Fig. 11, events are scheduled with 
respect to the master timing signal T^, (every 25 SM8's), the transfers AD and 
DA being initiated at the specified times TAD and TDA respectively. At the end 
of one computation cycle the digital program jumps to the main executive loop 
and waits there until T^, is received. In this manner, the digital section can be 
run synchronously with the analog. For instance, if T^ occurs every 50 ms and 
the time scale desired is 2. 0 times real time, the time step to be used in the 
digital program is J. sec. 
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2. Mode Control 

A very important function of the DOS 350 is to control the modes 
of operation of the system. Communication between the DOS and the 375 occurs 
through 16 sense lines which are set from the DOS and sensed on the digital 
section. In the other direction, eight flip-flops on the DOS can be set from the 
375 with special OCP instructions (output control pulses). Of these, four can 
also be reset from the digital console. 

All modes are controlled by pushbuttons from the DOS 350 and 
the following list summarizes the state and function of the analog and digital 
sections when the indicated pushbutton is depressed. 

(1) IC (initial condition) 

a. 231R in IC 

b. 375 in IC 

In IC the 375 goes through all computations with the 
exception of the integration routine. 

(2) TIC ( type in initial conditions) 

a. 2 31R in IC 

b. 375 is ready to accept new initialization data from 
the typewriter. From this mode the program returns 
automatically to the IC loop. 

(3) TYTI (type titles) 

a. 231R in IC 

b. The 375 types out title block and vline headings for the 
26 variables chosen for print out. 

(4) TRA (transfer only) 

a. 231R in IC 

b. 375 in IC 

This mode is for single stepping through the DA and AD 
transfers and was found very useful for problem checkout. 
The Master Timer is stopped and now incremented only 
manually by pushbutton action. The data is transferred 
in both directions, one word at a time, each time the 
pushbutton TEST CONV is depressed. After completion 
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of the DA transfer the digital program jumps to the 
executive waiting loop (EXWL) where it waits for another 
master timing signal from the appropriate pushbutton. 

DA Test Pattern 

a. 231R in IC 

b. The 375 goes directly to the DA transfer and back 

to the EXWL, skipping all calculations. A fixed block of 
data consisting of positive and negative maximum values 
(corresponding to ± lOOv on the analog) is transferred 
continuously. This was found very convenient for a quick 
check of the DA conversions. 

OP (operate) 

a. 231R in operate 

b. 375 in operate 

This is the normal mode of operation of the system 
TS (time scale) 

a. 231R in operate 

b. 375 in operate 

This pushbutton changes time scales, there being two 
arbitrary time scales available , i. e. , real time and 
twenty times real time. 

HOLD 

a. 231R in hold 

b. 375 in hold waiting loop 

The 375 stays in a waiting loop in the executive program 
until another mode is selected. 

DUMP 

a, 231R in hold during actual dump operation otherwise 
in IC or operate as previously selected. 

b. 375 goes to output routine at periodic intervals determined 
by TP the printout time interval. The digital jumps to the 
output routine and types out the present values of the 26 
variables. The DUMP command can be given in either IC, 
Hold or OP and the system resumes in which ever mode it 
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was at the time of execution. 

3. Data Transfers 

The data transfer for this problem is very demanding since two 
basic sampling rates are required, one for the short period aerodynamic 
functions and the other for the long period trajectory variables. The data 
transfers will be discussed in two parts: 1. ) the analog -to -digital conversion 
and 2. ) the digital -to -analog conversion. 

a. Analog -to -Digital Conversion 

Two variables, , and £ g are converted every 5 milli- 
seconds since they are used on the DOS for function generation. Once every 50 
milliseconds^, /3 , andyfcCare converted and transferred to the 375 for use in the 
long period trajectory calculations. Fig. 11 in the previous section shows these 
two different types of conversions. 

A block diagram of the ADC data transfer is shown in Fig. 13. 
Upon command from the master timing program on the DOS the ADC control 
logic increments the multiplexer to the proper channel and sends a convert 
signal to the AD converter. When the conversion is complete, the converted 
data is either loaded into a serial memory on the DOS for use in the function 
generation program or sent on to the 375 for use in the digital calculations. The 
data transfer from the DOS to the 375 is accomplished by loading the data into 
the buffer register and setting the parallel input channel ready flip flop on the 
digital computer. This flip-flop is enabled from the DOS by a pulse on the TEj^ 
line prior to each A-D transfer. The 375 then inputs the data through its parallel 
input channel. It is thus seen that the DOS controls all the A-D conversion, thus 
minimizing digital computing time on the 375. 

b. Digital -to -Analog Conversion 

Updated values of the aerodynamic force and moment 
coefficients which are generated on the DOS are transferred to the analog every 
5 milliseconds. The trajectory variables (altitude, temperature, guidance errors, 
etc. ) are transferred from the 375 to the analog once every 50 milliseconds. Fig. 11 
in the previous section shows the timing of these transfers. 

A block diagram of the DAC data transfer is shown in Fig. 12, The 
digital-to-analog transfers are initiated by the DOS control logic upon command 
from the master timing. If data is to be transferred from the 375 the DOS sets the 
output channel ready flip-flop on the 375. This flip-flop is enabled from the DOS 
by a pulse on the TE line prior to each D-A transfer. The 375 will then output 
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information to the buffer register and the DAC logic on the DOS loads it into the 
proper DA converter. To minimize digital computer time on the 375, four data 
words are loaded into the four DAC 1 s (under the control of the DOS). While this 
data is being converted and then demultiplexed on the analog under control of the 
DOS, the 375 is formating the next four words to be transferred. This process 
is repeated until all sixteen words are transferred. Using this technique, the 
total processing time consumed for the transfer operation (except for formating) 
is only 200yxsecs. The force and moment coefficients are transferred by loading 
the four DAC 1 s from the serial memories which store the latest computed values 
of these coefficients. The data is- then converted and demultiplexed on the analog 
and then the last four coefficients loaded. The DOS controls these conversions 
so that they occur at the end of each function generation period and therefore do 
not interfere with the conversions from the 37 5 to the analog. 

4. Function Generation 

/ 

/ 

/ 

The present simulation requires the generation of eight aero- 
dynamic force and moment coefficients. Among these are four functions of one 
variable j ^jt^( oL ) j CifsJ-*) ) O*<l(«0 ) > and four functions of two variables, 

are angle of attack, mach number, and elevator deflection respectively. These 
functions are generated on the DOS for the following reasons: 

(1) The functions of two variables are extremely difficult to 
generate in the analog section and would at best require a 
namber of sums and products of functions of one variable. 
The functions of one variable could be generated on the 
analog but were programmed on the DOS because of the 
ease of setup and the speed at which the functions could 
be changed to study other vehicle configurations. 

(2) The functions are also very difficult to generate on the 375 
because of the fast sampling rate required. The sampling 
rate on the functions should be at least 10 samples per second 
since they are used in the short period attitude loop of the 
vehicle. In the real time simulation the functions should 
therefore be sampled at least every 100 milliseconds and 

in the twenty times real time mode the functions should be sampled 
at least every 5 milliseconds. To meet this requirement 
and still use the 375, the 375 would have to be interrupted to 
compute these functions a great number of times during the 
major computation cycle. (In this problem the major compu- 
tation cycle is 50 milliseconds). It is estimated that the com- 
putation time required to compute these functions on the 375 
is about 4 milliseconds; therefore, if the program is stopped 
every 5 milliseconds for 
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function generation, only one millisecond of the five 
millisecond interval can be spent on the solution of the 
long period problem. Hence the major computation cycle 
would be 5 times longer and twenty times real time runs 
would be impossible. 

By use of the DOS, the function generation at sample rates of 
10 per sec per function can be accomplished in parallel with the digital calcu- 
lations in the 375. 

Fig. 14 shows a general block diagram of the function generation 
technique. The program for this problem utilizes 2-SM8's in series and allows 
32 curves of 16 points each to be stored. Linear interpolation between points is 
utilized, and for the functions of two variables several curves are used with linear 
interpolation between them. For the present problem the four functions of 
utilize 1 curve each and the remaining four functions of two variables are generated 
with 7 curves each. 

Each variable is sampled once per 2 SM8 cycle, i. e. , once every 
4 msec. Ld 20 times real time run this would correspond to a sampling rate 
higher thai 10 per second. 

5. Reaction Jet Control System 

Since portions of the flight are outside the atmosphere in regions 
where the dynamic pressure is too small to make use of aerodynamic surfaces 
for control, a reaction jet system is required. Fig. 's 15 and 16 show the block 
diagram of the digital controller, the jet configurations, the sign conventions and 
the combination of jets required to execute various commands. To make the 
reaction jet control easier, a pulse width modulation scheme is used which makes 
the moment proportional to the error signal. To prevent continuous pulsing of 
the jets even for small errors a deadzone is built into the controller. 

Fig. 's 17 and 18 show the circuitry required in the roll axis 
control to test the magnitude and sign of the roll error and to provide the palse 
width modulation . Fig. 19 shows the logic required to decode the commands so 
that the proper jets are used. 

The logic operations required in the simulation of the control system 
are performed on the DOS in parallel to all other operations in the digital section. 
The simulation of such a system by conventional analog techniques is a formidable 
task (dozens of switches and relays would be required). The use of the 375 for 
such an operation would result in a significant increase in digital computation time 
because of the large number of logic operations and the fast sampling rate required. 
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C. 



Digital Calculation on the 375 



The digital program was written to make the digital calculation 



time as small as possible and to stay withing a memory capacity of 4, 000 words. 
The material which follows gives a description of the equations which are solved 
on the 375 and the details of the digital computer programming. 



Fig. 20 shows the block diagram of the system of equations to be 
solved on the 375. A complete summary of the equations and symbols is given 
in Appendix A. 

The three degree -of -freedom tr anslational equations of motion are 
solved in a local horizontal coordinate systerr with axes along the north, east, and 
radial directions. (See Fig. 21) The gravity and altitude calculations are based 
upon an obla-e model of the earth, and the U.S. Standard Atmosphere, 1962, was 
stored in table form on the 375. Most of the equations are conventional and quite 
straighf orward, with the exception of the heat transfer, temperature, and guidance 
equations which will be discussed in more detail. 



1. Summary of the Digital Calculations 



The heat transfer and temperature equations are as follows: 
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is calculated from the stored density table by use of a second 
order curve fit technique to determine the derivative. These equations were 
calculated in the digital section because of the large numbers of complex functions 
required (i.e. , sine, cosine, log, exponential, square roots, etc.). However they 
are easily computed on the 375 by use of library subroutines. Fortunately these 
terms were also slow -varying, trajectory dependent quantities so that the slower 
serial computation on the digital is acceptable. The temperature rate equations 
involving angular rates are calculated on the analog because of thdr rapidly varying 
characteristic s. 



Because of their critical accuracy requirements, the guidance 
calculations are also performed on the digital section. These can be separated 
into two parts; (1) the determination of the range-to-go, (RTG) , heading angle 
to target, (CT^ ) , cross range-to-go, (CRTG), and heading error, and (2) the 
generation of down range and cross range errors (DR and CR^,) which are sent 
to the TRFCS and subsequently controlled to zero. The equations for range-to-go, 
heading angle, cross-range-to-go and heading error, which are determined from 
spherical trigonometric relations, are given in Appendix A. Fig. 22 illustrates 
the physical significance of these quantities. The down range and cross range 
errors used for guidance are generated as follows: 

DR E = RTG -RTG^ 
CR E = CRTG -CRTG A 
RTG A = f.(V) 



CRT<3A = -f z (V) 



f (V) and f (V) are the values of RTG and CRTG for a nominal re-entry. They 
are stored as a function of the relative velocity. 
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2. General Description of Digital Program 

In addition to solving the translational, temperature r and guidance 
equations, the digital program must accept mode control and timing commands 
from the DOS, scale variables as they are transferred in and out of the digital 
section, and finally, provide various input-output functions as described in the 
DOS section. 

Fig. 23 shows the flow diagram of the digital program. The 
executive program is discussed in detail in the section that follows. The dashed 
lines in the diagram link together the calculations which are performed when the IC 
mode of operation is selected. All operations are performed except the Runge-Kutta 
integration shown in block 9. When the operate mode is selected the integration 
loop is entered and three passes around the loop are made (see description of 
numerical integration which follows). At the end of three passes, the positions 
and velocities, temperatures and guidance functions are updated and the analog- 
to-digital and digital -to -analog transfers are made. The executive loop is then 
re-entered and the cycle continued until a different mode is selected. The following 
section discusses this executive loop in detail 

3. Digital Executive Program 

The main media of communication, through which the 375 receives 
the mode commands and timing from the DOS 350 is the digital executive program. 
(See Fig. 24). This program is essentially a chain t>f test instructions, through 
which the digital section interprets the DOS 350 commands and executes them by 
jumping to the respective portion of the stored program. 

Since some of the symbols are not identified in Fig. 24, a short 
explanation of the executive program follows. The first test determines whether T 
is greater than TF, where T is the current time (time elapsed sine the beginning of 
this simulation run •) and TF is final time indicating the end of the run. If final 
time is reached, results are outputed and the digital section is halted. "When re- 
started, the digital program begins with the next decision. The next decision 
determines whether fast or slow time scaling is requested. If fast time scale is 
selected, the increment size for integration DLT (delta t) , becomes 1 sec; otherwise 
it is 0. 1 sec. The subsequent two test instructions simply switch the digital 
program into TIC or TYTI modes if so commanded by the mode control logic. 
(See Section B for a definition of these modes). The T> TT decision determines 
if it is time for the periodic dump of pre -selected digital parameters. The variable 
TT holds the time for the next dump, say 20,30,40 etc. seconds. TP is the print- 
out time interval and increments TT when the dump time is reached. If the dump 
switch is on and T > TT the output dump is performed, if not, the program proceeds 
to the next decision. The next test, keeps the entire digital program under the 
timing control of the master timer. The 375 caanot proceed further until the next 
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Tl pulse indicates the beginning of the next computational cycle. Further testing 
takes place only after the Tl pulse arrives. The TDA test enables the programmer 
to place the digital program into a loop where only DA transfer is performed. This 
mode . serves as a convenient test for DA conversion. (This instruction may be 
removed from the program after the check-oat phase is over). The last few tests 
are self-evident and require no explanation with the exception of the seemingly 
superfluous second test for dump. The test for dump in hold mode makes it possible 
to dump any time by putting the computer in hold mode prior to the dump request, 
in addition to or in lieu of periodic dumps. In this manner the programmer can 
determine parameter values with digital accuracy at any time during the simulation. 

4. Numerical Integration 

Unquestionably one of the most complex and time-consuming 
parts of the digital program is the solution of the six simultaneous differential 
equations which provide the position (altitude, latitude, and longitude) and the 
velocity (radius rate^elocity east, velocity north) of the vehicle. 

The numerical techneque selected for this calculation was the 
fourth order Runge-Kutta method. The basic method as applied to a single differential 
equations is described briefly. 
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The increment for the second interval is computed in a similar manner by means 
of the same formulae 

A computer flow diagram for this scheme is shown in Fig. 2 5. The 

loop on the left side of this diagram corresponds to the constant calculations (k^, k^ 

k^, k^), while theloop on the right hand side of the drawing updates the independent 

variable and starts the calculation of the next increment. When the number of 

increments computed (n) is equal to the number desired (nf), the computation ends. 

Note that x and y , the starting values for the current increment calculation, must 
s s 

not be destroyed daring this process until the new starting values are produced, since 
all formulae are dependent upon them. 

The preceding technique was extended to a system of six equations 
in an obvious manner. 

While preparing the digital program and effort was made to combine 
the calculations due to the integration with other necessary calculations to minimize 
the processing time as well as the memory space requirements. This approach 
makes it difficult to trace the integration on the flow diagram of the digital program 
(Fig. 2 3). For example, the Runge-Kutta constants are determined in block 6 and 
the variables are incremented for the calculation of the next constant in block 9 
(such as Ve Ve ) e tc. , and the equations are evaluated with the incremented 

variables in blocks 3,4,6, and 7. The final calculation of the variables at 
is performed in block 9. 

The question may be asked by the competent reader: Why was 
th: s particular scheme chosen from many others available? In order to answer 
this question, the selection criteria Is listed below. 

(1) Self - starting method. 

It was desirable to choose a self - starting method to simplify 
programming and reduce memory requirement. (If the method 
is not self - starting, another method is needed to calculate the 
first few points of the solution. This virtually doubles the 
integration program). 

(2) Accuracy. 

The error of the fourth order Runge-Kutta method is of the 
order of h , and provides a wide enough range for h within 
which the accuracy of the calculation is acceptable. (During 
the simulation runs this reasoning proved to be a valid one). 
The accuracy of calculation using the second or third order 
method was considered to be insufficient or marginal. 




FIG. 25 BLOCK DIAGRAM OF NUMERICAL INTEGRATION 

(FOURTH ORDER RUNGE-KUTTA METHOD) 
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(3) Large time steps in integration. 

While processing time associated with this method is con- 
siderable, the accuracy is good enough to allow larger time 
steps. The increase in computation time using Runge-Kutta 
is more than off set by the increase in allowable step size. 

(4) Ease of changing step size . 

i In the selected method the length of the step can be modified 

at any time in the course of the computation without additional 
labor. This was considered to be a substantial advantage 
since the simulation is frequently called upon to change time 
scale, thus requiring a different time step for the integration. 

5. Utility Programs 

In order to satisfy the various data handling needs (type in, type 
out, punch tape, convert binary to decimal, etc. ) an excessive amount of digital 
programming must be done. Fortunately, all these programs are already available, 
tested, and clearly described in the 375 ! s "software package". This is no small 
feat if one considers that the above programs, together with the numerous sub- 
routines, normally add up to about 75% of all digital programming required. Further 
more, as it will be pointed out later, several programs are available to make the 
debugging and updating procedures efficient and fast. 

D. Analog Section 

The analog computer is the one link of the HYDAC 2400 system 
ideally suited for control system simulation by virtue of its capability for high 
speed, parallel computation and its input/output flexibility. Output data can be 
displayed in a multitude of forms such as X-Y plots, strip chart plots, oscilloscope 
displays, auxiliary meters, etc. Special purpose input equipment can easily be 
adapted for compatible operation with the analog computer. 

1. General Mechanization 

The block diagram of Fig. 26 delineates the mechanization of the 
analog program. A listing of the equations simulated in the analog section is 
given in Appendix A. The symbols used in the equations are defined in Appendix B. 

As previously mentioned the TRFCS and rotational dynamics were 
programmed on the analog due to their rapidly varying characteristics. When the 
vehicle is in the atmosphere, the vehicle attitude is controlled by aerodynamic 
control surfaces. The aerodynamic control moments are calculated from the surface 
deflections and the moment coefficients generated on the DOS 350. Out of the 
atmosphere, the vehicle attitude is controlled by a reaction jet control system. 
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The jet pulse logic is generated on the DOS 350 and transferred to the analog 
where the moments are produced and fed into the angular acceleration equations. 
The angular rates, generated from the angular acceleration equations, are used 
to calculate the Euler angles O and <f> . .These angles are used to resolve gravity 
into the body axis for use in the &t and (3 equations shown below. 

W rvi v ' £w 

W rvi 



These seemingly redundant force equations °^ and (3 are computed 
on the analog since they are required in the short period sensor equations shown 
below. The aerodynamic force coefficients in the OC and P equations are 
generated on the DOS 350. The geometric definitions of Qd and £3 are shown in 
Fig. 27. 

The short period sensor equations are: 



f = (\ -.l875o^) - .375 oCotT 



I is the temperature rate at the nose sensor; and/i"Tis the 
temperature rate differential between the two wing sensors (See Fig. 1). The 



V 



x<jz 



80Dy 



D£F INIT, 0N 



OF OC 



Hybrid Seminar 
Page 45 



stagnation point temperature information T and T which is used in these 
sensor equations is calculated in the digital section. The T and^AT are used 
in the TRFCS control equations shown below: 

S* = f,-(f ) 

S a = K, O-^c^ + X 2 P + K 3 £ R 

The T and T terms supply damping to the control equations 
alleviating the heating problems associated with an undamped trajectory. Closed 
loop guidance is achieved by adjusting the pitch axis controls with a compensated 
down range error, DR-rp » and adjusting the^^g to compensate for the cross range 
error, CR^,. The piloT can manually control the range by adjusting his temperature 
rate profile to eliminate the displayed down range and cross range errors. 

2. Cockpit Simulator and Display Equipment 

A cockpit simulator is utilized to evaluate the TRFCS in the 
manual mode, and is trunked directly to the analog computer. Computer outputs 
drive display meters on the TRFCS CONTROL PANEL* which monitor the following 
parameters: 



(1) 


temperatures 


(2) 


temperature rate 


(3) 


velocity 


(4) 


acceleration 


(5) 


pitch damper 


(6) 


yaw damper 


(7) 


bank angle 


(8) 


heading an.gle 


(9) 


elevator trim 


(10) 


range 


(11) 


range error 



* Refer to Fig. f s 28 and 29 of Reference #1. 
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(12) cross range 

(13) cross range error 

(14) altitude 

As seen in the analog program block diagram , many of these 
parameters are transferred from the digital section. Some of the parameters, 
though not necessary for control, indicate trajectory status and therefore main- 
tain the pilots confidence in his control information. The above parameters and 
other pertinent data are recorded on strip charts and X-Y plots for permanent 
record of each flight. 

Vehicle attitudes commanded by the flight control system are 
displayed on a large oscilloscope. An illustration of the display at various 
attitudes is shown in Fig. 28. 
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THE COMPUTATIONAL SPECTRUM OF HYBRID APPLICATIONS 



I. Introduction 

Several different approaches might be considered in classifying 
hybrid applications. The physicist would probably like to classify them ac- 
cording to basic physical law, such as conservation of momentum or the 
second law of thermodynamics. On the other hand, mathematicians would 
prefer to use as argument the various classes of equations, such as trans- 
cendental, or hyperbolic partial differential. These approaches adequately 
classify scientific "calculations" but would prove wholly impractical for that 
portion of scientific computation which is being considered here -- simula- 
tion. A given problem in systems engineering contains too many physical 
laws, too many types of mathematics for the above approaches. Therefore, 
this reduces to two other arguments -- the individual physical system (or 
application) and the computational nature of the simulation. The former 
would be useful, but would bear little relation to the types of computational 
equipment involved. Since the practicing engineer must give prime considera- 
tion to equipment available, the economics involved in the application compu- 
tational nature will be considered the most appropriate manner of classifying 
the spectrum of hybrid applications. 

II. The Computational Spectrum 

Now, consider the extent of the applications spectrum as shown 
in Fig. 1. Pure analog computation is on the left and pure digital on the right. 
Moving from left to right, the computational nature of the application will re- 
quire more and more digital techniques. The spectrum can be classified into 
eight distinct groups, each of which are associated with certain classes of ap- 
plications and equipment. These groups can be listed as: 

1. High Speed Sequential /Iterative Computation 

2. Data Analysis 

3. Function Generation 

4. Partial Differential Equations 

5. Optimization 

6. Digital Subsystems 

7. Dynamics /Extensive Algebra 

8. Dynamics /High Resolution 




FIG. I THE APPLICATIONS SPECTRUM 
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Group 1 corresponds to applications involving electronic mode 
control, point storage, one bit communication (comparators and DA switches) 
and a small amount of logic, as is found in the 231R-V. As these problems 
grow to involve many parameters or more elaborate measurement, the ex- 
tensive logic and one bit of communication of the DOS 350 becomes neces- 
sary, and the HYDAC 2000 becomes the appropriate computing system. 

Groups 2 through 5 require the addition of serial memory, whole 
word communication, and serial arithmetic as found in the more fully expanded 
DOS 350. These applications form the bulk of the requirement for an analog 
and a parallel digital section -- the HYDAC 2000. 

Finally, groups 6 through 8 require the additional capability of a 
serial digital section, and become the applications for the full HYDAC 2400 
system. As described earlier, the DOS 350 is the control center of the system. 
Therefore, the serial section is relieved to a large extent and can perform, 
at minimum iteration times, the calculation of algebraic equations and differen- 
tial equations at high resolution. 

A. High Speed Sequential /Iterative Computation 

The analog computer with a small amount of logic and point 
storage is in a position to automate many of the time consuming tasks of the 
past as well as open up new dimensions of system investigation. (See Fig. 2) 

1. Parameter surveys, the first step in non-linear optimization, 
can be completely automated with the use of electronic mode control, com- 
parators, D/A switches, track/ store accumulators, counters, and gates. 
Nested accumulation can provide automatic scan of many parameters, the slow- 
est changing under control of tape -set pots, while interpolated plots of variables 
vs. parameters, or vice versa, can be rapidly generated by use of strobing 
techniques. 

2. Monte Carlo methods are closely akin to the preceding, where 
highest speed accumulation can be a series of runs with noise inputs. Here the 
outputs are likely to be means or variances plotted against parameters. As an 
example of 1 and 2, refer to Fig. 3 . 

This is a block diagram of a high speed, two parameter survey 
of a system studied in the presence of noise. The system simulation might 



I. PARAMETER SURVEY 



L A 



2. MONTE CARLO 



— o 



n 



A- 






r 


I" 







X 



3. SPECTRUM ANALYSIS 




"e 



4. BOUNDARY VALUE 



5. OPTIMIZATION 



6. ALTERNATING CALCULATIONS 



FIG.2 HIGH SPEED SEQUENTIAL/ ITERATIVE COMPUTATION 



SLOW 



PARAMETH— H 
SCAN C 



MASTER 






SLAVE 




TIMER 




TIMER 





FAST 
PARAMETER 
SCAN 



SYSTEM 



SIMU- 



LATION 



EVALUAT- 
ION 



STROBE 



9i 











INTER- 




POLATION 





NOISE 



INTER- 
POLATION 



FIG. 3 MONTE CARLO PARAMETER SURVEY- BLOCK DIAGRAM 



Hybrid Seminar 
Page 6 



be a mechanization of a homing missile, where the integrators are under 
control of the Master Timer signals. One run might last ZO milliseconds. 
An ensemble of 50 runs is then counted by the Slave Timer, whose output 
(T s ) provides an advancing pulse to the Fast Parameter Scan accumulator. 
The accumulator, or ratchet circuit, is the basic circuit behind analog coun- 
ters, scanners, and evaluators, and is mechanized wi th a pair of track/ store 
amplifiers and a comparator. This Fast Parameter Scan circuit provides 
a series of levels, \\ , which can be forcing functions, initial conditions, or 
multiplier inputs in the simulation. The other output, T ss , is a digital 
signal for advancing the Slow Scanner. Thus, the \^ , for example, can be 
10 values of a loop gain, changing value every second, and the jj» 10 values 
of dead zone changing value every 10 seconds. Evaluation in this example 
might be measurement of the mean miss -distance over the 50 Monte Carlo 
runs, denoted as y* (t s ). This fast moving information must be strobed by 
a track/ store unit followed by a store /track unit. The output of the latter 
is a staircase, each level corresponding to the results of a given pair of 
the parameters X* L . Finally, an integrator and two amplifier circuit 

provides a straight line interpolation between levels of , when the 

abscissa, sayXi ^ s similarly treated. The result is a. family of curves, 
each curve representing response as a function of the parameterX^ for one 
value pi , and plotted out in about ten seconds. 

3. Spectrum Analysis is also a related topic, where now 
frequency is scanned with Nyquist or Describing Function plots as a result. 

4. Boundary Value problems offer a prime application for 
iterative techniques. Often, as with the tubular reactor, the relevant partial 
differential equation can be written as an ordinary differential equation with 
a final value constraint. In this example, the initial derivative is guessed 
and the final value measured. The steepest descents algorithm opei-ates on 
the resulting error and points the way to a better guess. Results of this 
series of iterations can become inputs to a real time control problem. 

5. Optimization over one or two parameters can be handled 
with this type of equipment, as can the related curve fitting problem. Errors 
are completed as the difference between output and a criteria function, or a 
least squares formulation, and operated on by the steepest descents algorithm. 
Manual and semi-automatic iterations are appropriate if there are more than 
two degrees of freedom. Here, point storage provides a visual display of the 
error, which is to be minimized by successive manual adjustments of para- 
meters. 
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6. Alternating Computation affords a convenient technique 
in certain cases, where the results of one computation will provide initial 
conditions for the next computation, and vice versa. 

B. Data Analysis 

The intent of this category is not to compete with well 
developed digital programs that now exist. Rather, it is the "quick look" 
aspects that are of interest, as so much of magnetic tape data analysis in- 
volves culling the good data from the bad. Also, there are a lot of simple 
tasks around that are wasteful of general purpose digital computer time. 
This category breaks into two areas: Simple tasks performed efficiently, 
and difficult tasks performed approximately, as indicated on Fig. 4. 

1. Logical Operations - defined as simple tasks of counting, 
measuring, formatting. Typical operations are: 

a. Counting number of times a magnetic tape signal 
crosses a threshold, or exists within given ranges, versus precise time 
marks. Comparators, flip-flops, gates and counters are the major tools. 

b. Peak reading within given intervals requires a small 
amount of serial memory, and a printer or punched paper tape in addition 
to A/D - D/A conversion. 

c. Formatting for IBM compatible tape, depending on rate 
requirements, can require a fair amount of logic, along with short, medium, 
and long delay lines for quick access bulk intermediate storage. 

2. Continuous Data Analysis - defined as analog techniques 

for obtaining statistical parameters and functions. The Exponentially -Mapped- 
Past approximation, and high speed control of integrators are the major tech- 
niques. Electronic mode control and logic are required for integration over 
precise intervals and control of magnetic tape. Possible parameters and 
functions are: 

a. Statistical parameters such as mean, variance, standard 
deviation, correlation coefficient. 

b. Statistical functions such as power spectral density, 
Fourier analysis, auto- and cross -correlation, LaGuerre function analysis. 




FIG.4 CONTINUOUS MEASURE OF POWER SPECTRAL DENSITY 
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c. Linear and non-linear regression. 

3. Communication control can vary from simple latching of 
comparators, through tape command for an A/D or D/A conversion, to free 
running A/D - D/A conversion. 

C. Function Generation 

In the past the analog computer has adequately handled 
most problems in function generation, through tapped potentiometer, diode, 
and integration techniques. But increasingly, simulations are requiring 
more extensive sets of functions with one, two and three variables, increased 
fidelity of transport delay, and more stringent resolution requirements, re- 
quiring digital techniques. Fig. 5 indicates the 3 basic categories. 

1. Empirical Functions, as required in two and three dimen- 
sional terrain avoidance, aerodynamic vehicles, etc. Serial memory, logic, 
serial arithmetic, and conversion provide the most efficient approach, par- 
ticularly when the physical system includes on-board radar. When a serial 
digital section is required for reasons such as algebra or resolution, function 
generation can be done in this unit when it is consistent with appropriate itera- 
tion times. It is not always appropriately handled in the serial digital section 
for the following reason: When the serial section becomes involved in generat- 
ing functions for the short period portions of the problem, total digital compu- 
tation for one iteration can become excessive. Here it is often appropriate to 
perform at least the short period function generation on the parallel digital 
section. Conversion control in its simplest form would normally permit a 
planned sequence of data transfer every iteration cycle. This sequence would 
have to communicate with both core and serial memory in the case of mixed 
function generation, thus versatility in the control circuits can be quite im- 
portant. The reader is referred to Chapter 3 for an example of function genera- 
tion using the parallel digital section. 

2. Transport Delay is conveniently and appropriately handled 
by serial memory, using the "marbles in a pipe" analogy. Time delay- 
frequency products of 30-60, and variable transport delay, are conveniently 
available. Conversion control would give a precession rate proportional to an 
analog variable, i. e. , flow. 

3. Cordic Logic can provide very fast generation of trigono- 
metric functions at a very fast rate using logic, serial memory, A/D-D/ A 
conversion, and serial arithmetic. One forward or reverse resolution may 
be performed every 500 microseconds, 2 resolvers may be implemented at a 
1 kc word rate, etc. Fig. 6 gives a diagram illustrating the basic definition 




FIG. 5 FUNCTION GENERATION 
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that is used; equation 1. Equations 2 and 3 follow from 1, and indicate the 
relationships used in the technique. 

a. Cartesian-To-Polar Conversion 

Sampled values of x and y are loaded into their 
respective registers, and k selected so that y is reduced. The number y 
can now be reduced to within one bit of zero in 16 steps, or 128 microseconds. 
This leaves a number proportional to R in the X register. The angle is ob- 
tained by accumulating values of _+<*i from delay line storage, determined by 
successive values of k. 

b. Polar-To-Cartesian Conversion 

The vector is rotated in successive steps by choos- 
ing the sign of k such that the original angle is reduced to zero in sixteen steps. 
The computation involved is essentially the inverse of the above casa 

D. Partial Differential Equations 

The analog all-parallel approach and the digital all -serial 
approach both show disadvantages in economical solution of PDE's. The all 
parallel approach requires too much equipment for solution for more than 10-20 
nodes. The approach generally used on the digital computer, the implicit solu- 
tion, requires solution of a large set of simultaneous algebraic equations for 
every time step, and is thus expensive. Various serial-parallel approaches 
show promise in providing accurate fine-mesh solutions in reasonably short 
times. 

1. Serial solution reverses the order of the all-parallel tech- 
nique. Time is discretized and the space variable made continuous, allowing 
simple iteration to meet a final boundary condition. (See Fig. 7). A function 
must be stored for the previous time step, implying the use of electronic mode 
control, logic A/D-D/ A conversion, and serial memory. The tradeoff now is 
time against mesh size. Conversion control must allow full rate sampling 
synchronized with the operate period. 

2. Hybrid Implicit Solution provides an "analog subroutine" for 
the lengthy matrix inversion required with the implicit numerical solution. 
The high speed arithmetic unit and core memory provides numerical integra- 
tion. Conversion control can provide a normal A/D-D/ A sequence every itera- 
tion cycle. 




FIG. 7 SERIAL SOLUTION OF DIFFUSION EQUATION 
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E. Optimization 

This topic may be considered to cover two allied areas - 
model development and adaptive control. The mathematical techniques are 
similar, making heavy use of steepest descents. 

1. Multiparameter Optimization becomes an extensive prob- 
lem in logic and control. In the general case of non-linear systems, coarse 
parameter surveys can first be made to assure an absolute and not a rela- 
tive extremum. Fig. 8 is a diagram of the steepest descent path and the 
quantized steepest descent path for a two parameter space. Mechanization 
of a quantized steepest descents formulation, assuming an appropriate 
criteria^ function, can be accomplished with electronic mode control, logic, 
comparators, and D/A switches. This can yield the speeds required to 
make optimization over many parameters a feasible matter, and a program 
versatile enough to readily adapt to any simulation. See Fig. 9 for a typical 
flow chart for an optimization program. 

2. Function optimization requires the addition of A/D-D/A 
conversion and serial memory to iterate on functions until certain criteria 
are met. This is similar to the solution of integral equations, where a 
kernel function is guessed and integrated iteratively until satisfying the equa- 
tion. 

3. Model development is becoming an increasingly important 
area of activity. Human factors, biomedical systems, hydrologic systems, 
econometric models - all are characterized by the need for synthesizing 
known output data into the appropriate mathematical models. Thus the known 
output becomes the criteria function, leading to the proper parameter set by 
use of the steepest descents technique. 

4. Adaptive Control again uses similar techniques, with the 
reference model as a method of generating criteria functions. This implies 
a continuous series of optimizations as external conditions change. 

F. Digital Subsystems 

One of the most obvious applications of hybrid computing is 
that of including digital devices as part of a dynamic model. Sometimes 
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thought of as sampled data control systems, they run the gamut from a 
small digital controller in a process unit operation to a completely integra- 
ted avionic system where all phases of operation, from communication and 
navigation to fire control and diagnostics, are under central control. 
Fig. 10 shows in block diagram form how such a simulation might be divided 
among the three basic types of equipment. Here the high speed arithmetic 
section makes its first appearance as a required element. At the same time, 
much care must be taken not to overload it with tasks more efficiently per- 
formed elsewhere. Therefore, radar simulation, function generation, and 
special purpose digital devices such as reaction jet control systems, are 
simulated in the parallel digital section, leaving the serial digital section 
free to simulate the central computer. Conversion control must be quite 
versatile to provide linkage to serial memory as well as to core memory, 
and put no restrictions on groupings or rates of transfer. In general, cer- 
tain channels must be sampled at maximum rate while others are sampled 
at various reduced rates. A significant amount of one bit communication is 
also necessary, such as comparators, D/A switches, sense lines, and out- 
put control lines. 

G. Dynamics /Extensive Algebra 

At times the major reason for choosing a hybrid approach 
lies in the impracticality of handling large sets of algebraic equations in 
dynamic systems. A case in point is the study of transient radiation effects 
in electronic circuits. Devices containing more than about two transistors 
have too many algebraic loops for practical all-analog simulation, yet the 
all-digital simulation is expensive. Retaining integrations on the analog sec- 
tion, algebra and control in the digital section, leads to economic simulation 
of very complex situations. Fig. 11 illustrates the division of computing 
tasks. Conversion control again should be versatile enough to handle certain 
channels at higher sampling rates than others. 

H. Dynamics /High Resolution 

This application may well be the classical one for hybrid 
computation. Simulation of space vehicles places such a contrasting demand 
on simulation equipment for speed and resolution that hybrid approaches are 
mandatory for economic solution. The concept here is to allow the serial 
digital section to provide the environment, i. e. , trajectory and dynamic pres- 
sure, as well as guidance, diagnostics, mission data, etc. The vehicle 
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attitude, control system, displays, etc. , are simulated on the analog section. 
This application in general overlaps the Digital Subsystems grouping in that 
digital controllers such as reaction jet control systems are generally present. 
The predominant need, however, is that of high resolution in the trajectory 
calculations for journeys to the moon and planets. Conversion control, as 
in previous sections, should provide no limitation as to sampling rates for 
individual channels, groupings of A/D and D/A transfers, and allow for a 
large amount of one bit communication such as comparators, D/A switches, 
sense lines, and output control lines. Chapter 3 gives ample illustration of 
this type of simulation. 

III. Summary 

As discussed above, the spectrum of hybrid applications is a 
broad one, and one that is developing quite rapidly. With software becoming 
available in each category, the development of these and new applications 
will show even faster growth. The user knows what applications are of prime in 
terest for his company now, and must govern himself accordingly. However, 
he usually does not know what might be vital to his company tomorrow. The 
user's main criteria should be to remain versatile and flexible, and use, to the 
greatest extent, the software and experience of the manufacturer and fellow 
users.,, and be aggressive in training personnel and educating management as 
to the benefits derived from hybrid simulation. 
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SOFTWARE IN HYBRID COMPUTATION 



I. Introduction 

Today when a new digital computer is announced and marketed, the 
manufacturer of the equipment is required to advertise and deliver, as part 
of his product, automatic programming language systems and various program- 
ming aids if he expects to be competitive. The original need for these program 
packets with written procedures describing their use developed when it was 
found at an early stage that the serial nature of these machines present inor- 
dinate difficulties for the applications engineer. Coding, program testing and 
interpreting large volumes of digital results proved most tedious, time- 
consuming and expensive. Some very clever techniques have since been em- 
ployed to reduce these undesirable activities. By using the computer to help 
in the preparation of programs (automatic programming), by taking advantage 
of previously written codes (library or canned routines), through automatic 
monitoring of a program during their execution (dumping, tracing, dynamic 
debugging) and by having the computer handle the needs of many users during 
a single time period (monitors, executive routines, an efficiency has been 
achieved through the creation of a programmed interface b etween the engineer 
and the machine. This interface, an adjunct to the hardware, is variously 
referred to as "software". This software development, therefore, by making 
digital computer programming easy for the non-computer -expert, is responsible 
for the almost universal acceptance of digital computer in the scientific field. 

The development of this sophisticated software has made it possible 
to increase computer utilization, to gain wider use of computers with minimum 
training of personnel, and to reduce duplication of programming effort for pro- 
grams of general utility. On the other hand, total dependance upon automatic 
programming, has the disadvantage of isolating the problem analyst, and even 
the programmer himself, from the computer. The analyst is restricted from 
communicating with his computer model while computation takes place. Simi- 
larly, the programmer is often limited in taking full advantage of the computer f s 
special features. 

In the analog computer field the reverse situation exists where no 
comparable "software" interface has grown up and the user communicates directly 
with the computer. There is little in the way of a requirement for the preserva- 
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tion of standard programs by analog computer programmers. As a result, 
the main feature of simulation by analog computation is that the problem 
analyst is directly involved in building the computer model and he main- 
tains rapport with his model from origination to completion of the problem. 

II. Hybrid Software 

The role that software must play in hybrid simulation of large 
complex systems is evident. "Hybrid software" must ease the programmer's 
burden, as it does for the digital programmer, and at the same time it must 
bring the analyst closer to his model rather than isolating him from it. 
Hybrid software must include not only coding for the sequential computer 
but also inter-connection diagrams and prewired patch panels for the paral- 
lel machine. The following types of software are needed to support growth 
of hybrid computation to meet the simulation needs of today. 

A. Compilers and Assemblers 

Conventional compilers and assemblers, providing automatic 
conversion from a symbolic to a machine language program, have a useful 
place in hybrid computation. In general, programming systems for hybrid 
computation may differ from conventional systems in four ways: 

(1) running time of the object program is minimized at the 
expense of compiling time. 

(2) actual running time for each program statement is 
precalculated or estimated to aid the programmer with 
timing aspects of the problem. 

(3) mostly machine -dependent programmer's language is 
utilized, since the programmer must be able to use 
special machine features to program for control of 
all interface operations. 

(4) the programmer's language is extended to all sections 
of the hybrid system. Such software, applied to an ef- 
ficiently designed hardware configuration, greatly 
facilitates the programmer's task in programming not 
only the digital, but also the analog and interface 
sections of a hybrid system. For example, setting 
analog function generators and potentiometers to pre- 
determined values, initializing integrators may be per- 
formed by some pseudo-instructions of the assembly 
language. 
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B. Debugging Aids and Diagnostics 

Conventional software aids for debugging a digital 
program include memory dump, updater and tracer routines, while diagnostics 
usually check out the proper operation of the digital processor hardware. 
Hybrid computer software, again, extends its capability to include debug- 
ging aids and diagnostics for the entire system. The following additional 
software of this type must be provided for efficient utilization of hybrid hard- 
ware: 

(1) Analog section diagnostics 

This is a symbolic program which accepts a problem- 
oriented language description of the analog flow 
diagram and produces automatically 

(a) setup documents (pot sheets) for the analog run 

(b) digital static check calculations for the analog 
program 

(c) possibly a dynamic check calculation for the 
analog program 

(2) Logic and Data Transfer Diagnostics 

This digital program is used to ascertain that the linkage , 
interface and logic hardware are operating properly. 
It checks out 

(a) the bidirectional data flow 

(b) the bidirectional control signal flow 

(c) the logic and secondary storage components, etc. 

C. Utility Library 

The conventional utility routines, prepared for a general- 
purpose hybrid computer include: 

(1) mathematical routines 

log, In, exponential, sine, cosine, tangent, arctangent, 
square root 

(2) computational routines 

single precision floating point, double precision floating 
point, double precision fixed point 
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(3) conversion routines 

(double and single precision) BCD/fixed point, BCD/ 
floating point, fixed point /floating point 

(4) input /output routines 

typewriter, paper tape, magnetic tape, automatic 
formatting on output. 

(5) algorithms 

a selection of numerical algorithms (integration 
methods, etc. ) 

(6) models 

subsystems, environments, etc. 

The last two categories require some additional explanation. 

In a hybrid system one may envision frequent need for various 
numerical operations such as definite integrals, solution of differential equa- 
tions, etc. Subroutines mechanizing various integration methods such as 
Runge-Kutta, Adams, Milne, etc. have been constructed. The different 
methods exhibit varying speeds and accuracies, and the user should choose 
to use those that fit best his problem. Furthermore, the hybrid simulation 
library should expand with routines for specific transfer functions of useful 
subsystem models that have become standard and are used in larger models; 
e.g., a typical servo controller. Another type of example is a function genera- 
tor program for any number of aerodynamic functions. A standard program 
for a complete aerodynamic vehicle simulation is- of value in a great number of 
aerospace computing laboratories. 

III. HYP AC 2400 Software 

EAI has had little need in the past to provide this type of support 
because of the parallel, engineer oriented aspects of the general-purpose analog 
computer. With the introduction of the HYDAC Z400 having, £.s a basic component, 
a straight digital computer EAI has chosen to improve its competitive position 
with respect to other manufacturers by: 

(1) offering a standard hybrid computer 

(2) providing comprehensive engineering backup and support 
for such a system 

(3) initiating development of hybrid- oriented automatic 
programming languages and aids. 



A. Basic Digital Software 
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1. Assembler - This is a metaprogram* which enables the 

user to write in symbolic machine language. A resultant translation provides 
the user with an edited listing of his original program and a condensed 
version of the code on punched paper tape. Several parts of a program may 
be "assembled" at different times later to be linked together and loaded 
automatic ally . 

a. DAP - (Digital Assembly Program) - a metaprogram 
based on the SHARE standard assembler. 



b. Macro-DAP - (for 8K-32K machines) - enables pro- 
grammers to define aggregates of codes under a single heading and thereby 
reduce programming effort. 

2. Interpreter - This is a monitor which executes pseudo- 
instructions to perform a given calculation. A complete new instruction re- 
pertoire is created for the engineering programmer which enables him to write 
in a DAP - like language but relieves him of the machine dependent aspects of 
DAP (input -output timing, scaling, storage allocation, etc. ) 

a. DIP - (Digital Interpretive Program) - a special 
translator -monitor which allows direct execution of symbolic program steps. 
Floating point calculation and input -output functions are automatic. 

3. Executive - a resident groups of subroutines which helps 
the programmer to test and correct his program dynamically. It will perform 
such functions as loading new programs, printing partial results, automatic- 
ally inserting or deleting instructions in the code, at the operator's direction. 

It may also contain standard input-output codes, interrupt bookkeeping routines, 
accounting routines and in fact all of those necessary for operation. 

a. DEP - (Digital Executive Program) 

(1) Relocating and Linking Loader 

(2) Input-Output Control 

(3) Editing and Monitoring functions 

(4) Debugging Aids (Dumps, Snaps, Tracer) 



* Metaprogram is a program used for the production of other programs. 
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4. Basic Subroutine Library - a set of subroutines which perform 

useful and standard calculations or manipulations vh ich can be easily incor- 
porated in the users own program. This library tends to grow as more and 
more diversified applications are programmed. A written procedure des- 
cribing the requirements, inputs, outputs, range of computation, limitations 
of each such subroutine accompanies a paper tape version (either in symbolic 
DAP, FORTRAN or machine code). 

a. Computational 

ADDX, SUBS, MPYX, DIVX - Double precision fixed point 
FADD, FSUB, FMPY, FDIV - Single and Double precision 
floating point 

POLY - Odd polynominal evaluation (used by many math 
subroutines) 



b. Mathematical 



(1) Fixed point - single and double precision 
SQD^ SQRX - Square Root 
EXPX, EXDX - Exponential 
LG2X, DPL2 - Logarithm (base 2) 
LGEX, DPLE - Logarithm (base e) 
LGAX, DPLT - Logarithm (base 10) 
SINX, SNDX - Sine X (input in radians) 
COSX, CODX - Cosine X (input in radians) 
TANX, TNDX - Tangent X (input in radians) 
ATNX, ATDX - Arctangent (output in radians) 



(2) Floating Point - single and double precision 



FSQR - Square Root 

FEXP - Exponential 

FLG2 - Logarithm (base 2) 

FLGE - Logarithm (base e) 

FLGA - Logarithm (base 10) 

FSIN - Sine (input in radians) 

FCOS - Cosine (input in radians) 

FTAN - Tangent (input in radians) 

FATN - Arctangent (output in radians) 



c. Conversion - single and double precision 



XBXD 
XDXB 
FBFD 
FDFB 



Binary -to -decimal (fixed point) 
Decimal to binary (fixed point) 
Binary to decimal (fixed to floating point) 
Decimal to binary (floating to fixed ;point) 
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FBXD - Binary to decimal (floating point) 
FDXB - Decimal to binary (floating point) 

d. Input -Output 

PTIN, PTOU - Paper tape I/O 
TYIN, TYOU - Typewrite r I /O 

Routines for on-line printer, magnetic tape, light pen/scope, disk are also 
available. 

e. Utility and Debugging 

DUMP - output to typewriter, paper tape, printer magnetic 
tape in octal, decimal, alphanumeric or mnemonic instruction format. 

DEBUG - Program testing communications package for on- 
line typewriter debugging - incorporates ability to modify or output data, search 
or clear memory, input correction tapes or output self loading dump tapes, set 
breakpoints, etc. 

TRACER - Enables an automatic step-by- step monitored 
execution of a program. Intermediate register results are printed on the type- 
writer or on-line printer. 

UPDATE - Enables the updating and correction of symbolic 
program tapes and/or data tapes. 

UPSUB - Enables updating of the tape containing the library 
of subroutines - inserts, deletes and duplicates. 

LOADER - Used to load absolute and relocatable DAP and/or 
FORTRAN object (binary output) taped programs. Automatically scans the 
library tape for standard subroutines and loads and links them to main program. 
It is the 375 equivalent of the SHARE 704/709-90 BSS Loader. 

MISC - Quick-Look-and-Change, Insert (formatting editor), 
Relocatable Dump/Punch (octal output), Memory Clear, Bioctal Dump for self 
loading tapes, Tape-to-Tape conversion (for tapes with a code different from 
375 input BCD code). One wo uld have to see the writeups for the above in order 
to gain insight into their use. 
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5. Compiler - a metaprogram which translates directly into 
machine code the statement of a problem. An artificial language created 
specifically for expressing easily certain types of computations and/or 
logical manipulations (i. e. algebraic or engineering problems), together 
with a language translator is implemented on a computer. Many problem- 
oriented languages (POL's) exist. One such language, FORTRAN (FORmula 
TRANslator) is accepted in the U.S. as standard and an implementation 
exists for almost every machine. (N.B. Unfortunately no minimum standard 
definition has been accepted for the FORTRAN language. All implementa- 
tions of it vary in scope for different computers, hence it is not machine- 
independent. Unless the programmer takes extreme care in using compatible 
subsets of the language(s) he cannot achieve equivalent results for parallel 
computations on dissimilar machines.) 

a. FORTRAN - An engineering/algebraic compiler, easy to 
use and learn, flexible automatic formatting of input-output, easy to debug. 
Object programs compatible with DAP condensed output tape format. 

(1) In- Line Machine Language Coding with FORTRAN. 

The 375 FORTRAN Compiler provides special provision 
for real-time and hardware -oriented programming. A capability to write 
DAP-like entries "in-line" with the algebraic statements gives the programmer 
a control over computation not afforded by any competitive FORTRAN system. 
This ability is a necessity if any compiler metaprogram is to be used to 
generate the codes for a hybrid c omputation. 

(2) Expanded features of 31 S FORTRAN. 

Although compatible with a generally accepted FORTRAN 
II, this 375 compiler has been extended to include the most desirable features 
of a new FORTRAN IV version as defined by IBM for its large scale computers. 

B. Expansion of 37 5 Software for Hybrid Computation 

1. Mathematical Library 

EAI has provided many diversified subroutines for 
types of numerical integration . (Runge-Kutta, Adams, Milne, Euler, Predictor- 
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Corrector, etc) which exhibit varying speeds and accuracies. The user 
may choose to use those which fit best his problem or he may incorporate 
several at once to be optionally selected during the course of a run. This 
selection may be under programmed or manual control. Routines for high 
speed function generation and non-linear transfer functions are also avail- 
able. 

2. Automatic Programming Aids 

HYTRAN (HYbrid TRANslator) will be used to refer to a 
family of metaprograms being implemented on the HYDAC 2400 computer to 
provide digital assistance in the programming of analog consoles. The three 
HYTRAN programs being written are: 

(1) An Interpretive Static Check Generator 

(2) As an extention of (1), a Documenting program 

(3) An On-Line Diagnostic Generator 

The first program reads and interprets the physical 
problem and the analog program description which is punched on paper tape 
in an analog -oriented (HYTRAN) language, computers parameters, initial 
conditions, and both the physical and the voltage Static Check. It examines 
the two Static Check values for consistency, pin-pointing any discrepancy or 
overload at the component level, and finally outputs an alphabetic list of static 
check values. 

The second generates complete documentation of the 
analog program, featuring Pot-Sheets, Amplifier -Sheets, Alphabetic lists of 
all Parameters and Variables, Cross -Reference, and Pot-Setting and a Static 
Check tape in ADIOS*format. 

The third and final program checks the measured Static 
Check Voltages for consistency with the analog circuit diagram and specifica- 
tions of the analog components, providing a rapid means of locating mispatching 
or component failures. 

Many additional HYTRAN metaprograms are planned in- 
cluding other variations on the APACHE, DAS and MIDAS efforts. 



* Analog Digital Input-Output System 
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3. HYDAC - ALGOL - a self -compiling compiler to be used 

for the definition of a command and control language for HYDAC Z400. This 
special type of translator has flexibility enough to enable the user to generate 
his own software with relative ease. Future expansions in the hybrid hard- 
ware can be immediately reflected in an updated version of this real-time 
programming tool. ( N. B . the implementation is to be based on the NELIAC 
system as developed for real-time command and control at the Naval Electronics 
Laboratory, San Diego, California) 

C. HYDAC 2400 Software - (Operating on 375/DOS/231-R-V) 

1. Dual-Processing Executive Program 

For hybrid systems with at least 8K 375, a provision is 
made through hardware modification for independent computation of two problems 
in tandem. The 375 memory and I-O equipment may be shared and a hybrid 
computation can continue at the same time a purely digital computation (i. e. a 
FORTRAN compilation) is being processed. No interference between the two 
co-residents will occur because the storage lockout /protect mode is controlled 
externally through the logic of the DOS. The Executive Program resides in 
375 memory and services the interrupts and performs other important house- 
keeping functions. 

This feature is very important where the checkout and exe- 
cution of a hybrid computation dominates the time available for using the 375 
for program preparation. By time -sharing on HYDAC 2400 a greater efficiency 
is realized over any competitive system. 

2. Modifications to DAP and FORTRAN for control of DOS/ 37 5 
interface and further control of 231-R-V are anticipated. These additions allow 
symbolic programming to include: 

a. Ability to perform ADIOS functions 

b. . Real-time iterative 231-R-V mode control through DOS 

c. Serial Data Transfer Unit* and basic system control 

d. 315 conversion system* control 

3. HYDAC Monitor System (for dynamic program debugging). A 
set of hardware and software procedures are developed for automatic program 
checkout. These include such debugging packages as: 



* Sections of the DOS 350 Computer 
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a. Tracer for 375 

b. DOS memory dump and selective printout 

c. AERO scan and 231-R-V readout and selective print 
, d. Automatic pot setting and checking 

e. GPAC Check Program execution. (This monitor will be 
incorporated with the DEP (Digital Executive Program) 
for standardization of on-line debugging techniques. ) 

4. Miscellaneous Hybrid Utility Programs - The se are being 
developed as their need arises. Some are already available. 

a. Tape-to-Tape Conversion (ADIOS/DOS/375) 

b. Digital Plotting Subroutines. Strip Chart, X-Y, CRT 

c. DOS/ 37 5 direct memory linkage software 

d. DOS as satellite digital computer for 375 for use as: 

(1) Digital Resolver 

(2) Auxiliary input-output device 

(3) Additional memory and logic capability for parallel 
computation 

(4) An ultra high-speed function generator 



e. Use of DOS as real-time interrupt and control computer 
for 375 during complex hybrid simulations 



IV Conclusions 



EAI, through its own growing experience in hybrid computation, 
provides the HYDAC users with the necessary programming tools. Software 
never has solved a problem by itself. It can be used to get the solution sooner, 
however. Unique facilities are provided through the metaprograms for the 
user to produce his own software -- often reflecting his applications and 
parochial needs. There is no end to this process. As knowledge in computa- 
tion grows, more difficult and even impossible, problem solutions are sought 
and solved. Software is improved, expanded and newly created in the cus- 
tomer's own laboratory. The preponderance of useful programming aids has 
always originated with users of the machines. This manufacturer hopes never 
to reverse this trend. By providing a comprehensive software package in the 
beginning, the hope is to get the users "on board" more quickly and to turn 
out useful work in short order. EAI points with ever increasing pride to the 
solutions which can be and have been implemented on the hybrid computer. 
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APPENDIX A 

I. Digital Comput er Equa tions 

A. Basic Translational Equations in Local Horizontal System 

The translational equations are solved in a local horizontal co- 
ordinate system with axes along the north, east, and radial directions. Fig. 22 , 
Chapter 3 show this coordinate system. 

V E ^fB (f sin.fl)V w ~(f COsQ)r (3) 

B. Aerodynamic Force Calculation and Transformation 

The aerodynamic forces are calculated with respect to the relative 
velocity vector {lift, drag, and y axis acceleration) and then transformed into the 
local horizontal coordinate system. 



D - <JSC D 



(4) 



(5) 



Y= ^5C Y (6) 



(7) 
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r 



+ .2. 1 CObO" (9) 



C. Angular Rate and Position Calculation 

The derivatives of the space fixed coordinates Q. and fare solved 
as functions of the inertial velocities V and V^.. The derivative of the longitude, 
is calculated by subtracting the earth's rotational rate from f . Since the earth 
rotates in the logitudinal direction CL is the latitude. Altitude is calculated with 
respect to an oblate earth. 

(10) 



rcosA 



(11) 



A = f - Wo (12) 



\n = r - r 6 + si^il d 3 ) 

2 0- e 2 ) 



D. Velocity, Velocity Heading and Flight Path Angle Calculations 

The total relative (to the earth) velocity is calculated from the 
velocity components in the N, E and radial directions. The horizontal relative 
velocity has two components, east and north. The easterly relative velocity is 
the difference between V and the tangential rate of the earths rotation. The 
relative northerly velocity is the same as the northerly inertial velocity. The 
flight path angle ^ , is the angle between the horizontal and the total velocity 
vectors. The velocity heading angle C, is measured from east to north. 
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V = V V e*+ V »W*' (14) 



\r + (is) 



sin * = ^ ( 17 > 



COS)! = ^ (18) 



if ~ tan* ~ 



V 



M 



(19) 



5m(T-- ^ (20) 



eescr-- -Ve* (2i) 



<T - tan 1 2b± (22) 



E. Density and Dynamic Pressure 



The density is tabulated in square root form to compress the 
range of variation and was taken from U. S. Standard Atomophere, 1962. The 
dynamic pressure is a function of density and velocity. 
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ip =^fp (K) (23) 



F. Gravity Calculation 

The gravity components are calculated for an oblate spheroid 
shaped earth. Since the earth is assumed to be rotationally symetrical, 

g E =o 



q u = -i^JlISi s 'm^aacssi (26) 



G. Heat Transfer and Temperature Equations 

These equations are generated in the digital program and trans 
ferred to the analog for calculating parameters in the flight control equations. 



i c - cos LS aC (29) 



(30) 
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1 

T, = I — + (3D 



s 



T = -* fa. + \ \ \ H (32) 



(34) 



>e **" -2T^T- (35) 



H. Target Heading Angle, Range-to-go and Cross Range -to-go 

The range and heading angles are calculated using spherical 
trigonometry and the resulting angles are multiplied by the earth radius yielding 
the range. 



T^A Q% - -r . (3b) 

SWA iA^ A) 



^ * e Krr= F 3m(AT-A) {37) 

"taiiCR = ~ ccr T - cr) ^ 6 KTC> , (33) 
CRTG = (CR) r n 
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I. Heading Error 

The angle between the great circle heading to the target, CT T , and 
the velocity heading, G~ . 



He.adittcj Error - O^-CT (39) 

J. Guidance Equations 

Error equations for use in the TRFCS 

DR E = RTG - RTG a (40) 
RTG d = -f.(V') 

CR e =CRT6-CRT6a (4i) 
CRTS d = i(V) 

II. An alog Comp uter Equa tions 
A. Angular Acceleration 

The following equations are calculated in the body system: 

P - -(l£Z±l\Q& + L * M « (42) 



V I* i 1 
Q= (l M - k\gp 4 M i My (43) 

j? = - / lllll 1 PQ ■+ N ± (44) 

V U J I* 
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B. Aerodynamic Torqj.es 

These parameters are generated by the moment of aerodynamic 



forces. 



L = <J S b [C /(S (oO 6 + C A 6 r * C^6«~] (45) 

H= 6.)] (46) 

N = <J S b[c„^)/3 + C^^r+C^w.)^ (47) 



C. Euler Angles 

Generated for local horizontal to body axis transformation. 



O = y 1 (Qcostf - Rs\n<t>)At 



(48) 



(49) 



D. Attitude Rate Equations 

Generated in this fjprm for use in the heat sensor equations. 

VV Vr / mv 

+ Q" ( P^^oc -ti?siN)i3 (50) 

V V Vr / v 7 tnv 

@ - ^3r" ^(cw^cw^/^-s/^^^jS (51) 
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^Cf ~ PCOS^ 4- G6+ R> S\V\ (52) 



E. Heat Sensor Equations 

These short period equations are used for control and are 
discussed in the text. 



T = t s (i- -I875<X 2 ) - .375o6*T;s (53) 



AT - # o2.i T s jS (54) 



AT - .02i (T s ^ -t- T 5 /3) (55) 



F. TRFCS Control Equations 

The control parameters include the variable surfaces and 
the commanded bank angle. 



<S B = -f.CT) (56) 



6 a = K. 0*—^ G ) + if^P + |C 3< £ R (57) 



6 R = -f z (AT) + (58) 



= •F 3 (T) (59) 
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APPENDIX B 
DEFINITION OF SYMBOLS 

Definitio n 

acceleration north (local level coordinates) 

acceleration east (local level coordinates) 
acceleration up (local level coordinates) 

aerodynamic acceleration north 

aerodynamic acceleration east 

aerodynamic acceleration up 
gravitational acceleration up 
gravitational acceleration north 

lift acceleration 

drag acceleration 

side force acceleration 

velocity north, inertial 
velocity east, inertial 
radius rate 

t otal inertial velocit y 
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D efini tion 

total relative velocity- 
relative horizontal velocity- 
velocity north, relative (equals V ) 
velocity east, relative 

mach number 

latitude rate 

longitude rate, inertial 

longitude rate, relative 

rate of change of angle of attack 

rate of change of angle of sideslip 

roll angular acceleration 

pitch angular acceleration 

yaw angular acceleration 

roll rate about the body x-axis 

pitch rate about the body y-axis 

yaw rate about the body z-axis 

Euler rates of body axis with respect 
to local horizontal 

aerodynamic torque about the body x axis 
aerodynamic torque about the body y axis 
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Va riable Definition 

N aerodynamic torque about the body z axis 

M reaction jet moment about the body x axis 

reaction jet moment about the body y axis 

M reaction jet moment about the body z axis 

latitude 

^ longitude 

^ angle of attack 

^ angle of sideslip 

bank angle (lift vector with vertical plane) 
$ body roll angle 

® body pitch angle 

¥ body heading angle 

bank angular rate 

elevator deflection angle 
6cl aileron deflection angle 

6 j~ rudder deflection angle 

^ flight path angle ( +■ for r-f) 

®~ velocity heading angle ( to the north from east) 

Or heading to target (+ to the north from east) 

air density 

V" radial distance from earth center to vehicle 
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Variable 



Definition 



h 

DR E 

CR E 
RTG 
CRTG 

C Y 

q 

c. 



M 

c, 



a. 



rs 



cs 



altitude 

down range error 

cross range error 

range -to-go 
cross -range -to -go 
side force coefficient 
dynamic pressure 
lift coefficient 

drag coefficient 

pitching moment coefficient 

rolling moment coefficient due to sideslip 

yawing moment coefficient due to sideslip 

rolling moment coefficient due to aileron deflection 

yawing moment coefficient due to aileron deflection 

yawing moment coefficient due to rudder deflection 

rolling moment coefficient due to rudder deflection 

radiative heat transfer rate at stagnation point 

convective heat transfer rate at stagnation point 

radiative heat transfer rate at sensor location 



convective heat transfer rate at .sensor location 
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Definition 

total integrated heat flux at stagnation 

temperature rate at sensor location 
temperature rate at stagnation point 

temperature at sensor 
temperature at stagnation point 

differential temperature rate between 
wing sensors 

earth gravitational constant 
earth oblateness factor 
radius of the earth -equatorial 

eccentricity of the earth 
rotation rate of the earth 
emissivity 

Stefan Boltzman constant 
sea level gravity 
mass of the vehicle 
wing area 
wing span 
wing chord 

moment of inertial about x axis 
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Definition 

moment of inertia about y axis 
moment of inertia about z axis 
nose cap radius 

exponential index of the atmosphere 



